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Abstract: Significant new opportunities for astrophysics and cosmology have been identified at low 

radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere 
designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with ar- 
cmiimtc angular resolution and high survey efficiency. The telescope will enable new advances along 
four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisa- 
tion in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain 
astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Wide- 
field Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) 
low-band telescope and is the only low-frequency ^KA precursor facility. In this paper, we review the 
performance properties of the Murchison Widefield Array and describe its primary scientific objectives. 

Keywords: cosmology: dark ages, reionization, first stars — instrumentation: interferometers — radio 
continuum: general — radio lines: general — sun: general 
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1 Introduction 

The Murchison Widefield Array (MWAQ is a low- 
frequency radio telescope operating between 80 and 
300 MHz. It is located at the Murchison Radio-astronomy 
Observatory (MRO) in Western Australia, the planned 
site of the future Square Kilometre Array (SKA[^ low- 
band telescope, and is one of three telescopes desig- 
nated as a Precursor for the SKA. The MWA has been 
developed by an international collaboration, including 
partners from Australia, India, New Zealand, and the 
United States. It has begun operations in 2013. In 
this paper, we review the science capabilities of the 
telescope and the primary science objectives motivat- 
ing its design and operation. Details of the technical 
design and specifications of the MWA as built are pre - 
sented in a companion article (Tingay et al. 20131, 
which supercedes the earlier Lonsdale et al. ( 2009 1 de- 



Table 1: System parameters for MWA 



scription of plans for the instrument. 

The four key science themes driving the MWA are: 
1) searching for redshifted 21 cm emission from the 
epoch of reionisation (EoR) in the early Universe; 2) 
Galactic and extragalactic surveys; 3) time-domain as- 
trophysics; and 4) solar, heliospheric, and ionospheric 
science and space weather. We describe the plans and 
objectives for each of the key science themes in Sec- 
tions 2 through 5. 

The frequency range targeted by the MWA is one 
of the more difficult portions of the spectrum for ground- 
based radio telescopes due to strong distortions of in- 
coming signals by the Earth's ionosphere and also due 
to congested terrestrial communications, including FM 
radio, TV stations, aircraft navigation, and satellite 
communications. Hence, the low-frequency radio sky 
remains one of the least explored domains in astron- 
omy. The MRO site mitigates many issues associated 
with radio interference for the telescope, but process- 
ing of MWA observations requires particular care to 
correct for the ionospheric distortions. 

The unifying design theme for the MWA is that 
of high survey efficiency. Many of the science goals re- 
quire large areas of the sky to be surveyed and the EoR 
objective further requires very high surface brightness 
sensitivity. It was recognized early in the development 
of the MWA that arcminute angular resolution is a 
convenient balance between the opposing tensions of 
maximizing scientific returns and minimizing techni- 
cal complexity. Arcminute angular scales are suffi- 
cient to probe typical ionised bubbles during the EoR, 
study diffuse and polarized emission structures in the 
Galaxy, reduce extragalactic source confusion, image 
bursts in the solar corona near the Sun's surface, and 
localize transient events with sufficient precision to al- 
low follow-on searches at other wavelengths. In order 
to achieve arcminute resolution, the longest baseline of 
the MWA is 3 km. This distance is within the scale size 
of a typical isoplanatic patch in the ionosphere, allow- 
ing simplying assumptions to be employed for widefield 
ionospheric calibration and keeping real-time imaging 
computations within reach of modest computer clus- 
ters. 



Parameter 



''Value 



Number of tiles 


128 


Tile area (m^) 


21.5 


Total collecting area (m^) 


2752 


Receiver temperature (K) 


50 


Typical sky temperature (K) 


350'' 


Field of view (deg^) 


610 


Angular resolution (arcmin) 


2 


Frequency range (MHz) 


80-300 


Instantaneous bandwidth (MHz) 


30.72 


Spectral resolution (MHz) 


0.04 


Temporal resolution (s) 


0.5 


Minimum baseline (m) 


7.7 


Maximum baseline (m) 


2864 


Estimated confusion limit (mJy) 


10 



"For frequency-dependent parameters, listed 

values are given at 150 MHz 
''Sky temperature varies considerably with 
G alactic latitude . Here , we use typical values 
from Nijboer et al. (2009) and Rogers & Bowman 
020081) 



The system properties of the MWA are described 
in detail by Tingay et al. (20131. We summarize them 
here in Table [Ij The time needed to reach a point 
source sensitivity of for the MWA is given by: 



2kBT 



1 

a?Bn„ 



(1) 



where fee is the Boltzmann constant, T = Tsky + Trcv 
is the system temperature, AcH is the effective area of 
each antenna tile, is the number of antenna tiles, 
tc is the correlator efficiency (assumed to be unity), B 
is the instantaneous bandwidth, and rip is the number 
of polarizations. For continuum imaging at 150 MHz, 
this reduces to i « 8 x 10'^ /{Baf) seconds for B in 
units of MHz and as in units of mJy. The intrinsic 
source confusion limit for the MWA angular resol ution 
of 2' is estimated to be ~ 10 mJy at 150 MHz ( Con- 



donl[T974l [Windhorst et aI][T985| [Condon et aLlpOSf 



^ http: / / mwatelescope.org 
•^http: / / skatelescope.org 



Subrahmanyan fc Ekers||2002| ). Hence, the MWA will 
reach the confusion limit in ~ 25 seconds for contin- 
uum imaging or within several hours in a single 40 kHz 
spectral channel. 

The surface brightness sensitivity of the MWA varies 
significantly with angular scale. The MWA antenna 
layout emphasizes a dense, central core in order to 
achieve a high surface brightness sensitivity on degree- 
scales for EoR science. For a typical EoR spectral res- 
olution of 1 MHz, the brightness temperature uncer- 
tainty is less than 150 mK after 1 hour of integration 
on degree-scales, whereas it is approximately 3 K on 
arcminute-scales . 

In the remainder of the paper, we discuss the four 
key science themes. Section 2 describes the epoch of 
reionisation science. In Section 3, broad categories 
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of radio astronomy science are discussed within the 
context of sky surveys. Section 4 describes the time- 
domain astronomy that will be possible with the MWA 
and Section 5 presents the primary solar, heliospheric, 
and ionospheric science goals. We conclude in Sec- 
tion |6l 



2 Epoch of Reionisation 

One of the most significant events in the history of 
the intergalactic medium (IGM) was the reionisation 
of hydrogen. The epoch of reionisation marks the time 
when the small fraction of matter within galaxies sig- 
nificantly affected their surroundings (and hence future 
generations of galaxies), ionizing the intergalactic gas 
and heating it from 10s of Kelvin to ~ 10* K. 

The MWA has been designed with the goal of de- 
tecting 21 cm emission from the reionisation epoch. It 
is one of three new radio arrays specifically targeting 
this science objective, along with PAPEErl and LO- 
The 21 cm hyperfine transition line of neutral 
hydrogen is anticipated to provide an ideal probe of 
reionisation. Its weak oscillator strength (in contrast 
to the hy-a line) allows 21 cm photons to propagate 
from fully neutral environments at high redshifts with- 
out appreciable absorption. This property makes it 
possible to directly observe the neutral hydrogen gas 
before and during reionisation. 

Unlike the continuum emission of the cosmic mi- 
crowave background (CMB), 21 cm emission contains 
information in the frequency (or redshift) dimension in 
addition to in angle, enabling the study of the full his- 
tory of reionisation. Redshifted 21 cm emission from 
neutral hydrogen gas in the IGM should appear as 
a faint, diffuse background in radio frequencies be- 
low v < 200 MHz for redshifts above 2 > 6, where 
ly = 1420/(1+2:) MHz. At these redshifts, the expected 
unpolarized differential brightness temperature, 5T21, 
of the redshifted 21 cm line relative to the CMB is 
readily calculable from basic principles (see 



(Barkana & Loeb 20051. The other two factors de- 



Furlanetto 

et al.|2006 for a complete review) and is given by 



<5T2i( 



1 - 



25 x„i (1 + 5) 
Ts 



1 + 2 
10 

H{z)/{i + zy 



di^l|/cir|| 



1/2 



mK, 



(2) 



where xhi is the neutral fraction of hydrogen, S is 
the local over-density in baryons, 2^7(2) = 2.73 (1 + 
2) K is the temperature of CMB at the redshift of 
interest, Ts is the spin temperature that describes the 
relative population of the ground and excited states of 
the 21 cm transition, H{z) is the Hubble constant, and 
di/||/dr|| is the velocity gradient along the line of sight. 

All of these factors contain unique astrophysical 
information. The dependence on S traces the develop- 
ment of the cosmic web | |Scott fc Rees|1990[ ) , while the 
velocity factor creates line-of-sight distortions that can 
be used to separate different components of the signal 



pend on the ambient radiation fields in the early Uni- 
verse. The local neutral fraction, xhi, is determined 
by the ionizing radiation emitted by galaxies and pos- 
sibly quasars. Theoretical studies suggest that the dis- 
tribution of neutral gas during reionisation is very in- 
homogeneous (e.g. McQuinn et al.|2007 l. Finally, the 
21 cm spin temperature, Ts, is determined primar- 
ily by three mechanisms that affect the hyperfine level 
populations: the ambient radiation at 21 cm due to the 
CMB with T-y, the ultraviolet background that tends 



to drive Ts to the gas temperature ( Wouthuysen|1952 



Field|[T958 ), and heating from the cosmological X-ray 
background (Che n fc Miralda-Esc udc 2004). By the 
time the Universe is a few percent ionized, most cal- 
culations predict that Ts ^ T-y, such that fiuctuations 
in Ts are not important (e.g. Furlanetto et al.|2006[ ). 

Figure [1] illustrates how the factors in Equation [2] 
combine to create the observable 21 cm signal. The 
panels in this figure are from the simulations of |Mc-| 
'Quinn et al. (2007) and are 100 cMpc across or ~ 0.5° 
in angle (a small fraction of the ~ 20° MWA field of 
view) . The left panel of Figure [l] shows the expected 
signal just before reionisation begins, when the fluc- 
tuations are dominated by local perturbations in the 
matter density. In the middle panel, reionisation is 
underway (xhi ~ 0.5), and the large-scale voids in 
the map signify ionized regions. After reionisation has 
largely finished (right panel), only localized pockets of 
neutral hydrogen remain. 

The MWA will attempt to probe structures in the 
redshifted 21 cm background in order to reveal the as- 
trophysics of reionisation. The MWA is expected to 
be sensitive to the power spectrum of the 21 crnsig- 
nal over the redshift range 6 < 2 < 10 (Section 2.1 1. 



^http:/ /eor. berkeley.edu 
■* htt p : // www . lof ar . org 



Its observations have the potential to characterise the 
nature of the sources that are responsible for ioniz- 
ing the IGM, chart the evolution of the global neutral 
fraction, and probe the nature of quasar emissions by 
constraining the properties of their ionized proximity 
zones. 



2.1 21 cm Power Spectrum 

Initial MWA observations aimed at detecting 21 cm 
emission will focus on two target regions at relatively 
high Galactic latitudes (see Figure [2]|. Over the course 
of one year, approximately 1000 hours of optimal night- 
time observations will be divided between these two 
regions. The signal-to- noise ratio (S/N) in each pixel 
will be < 1, with the highest S/N on the largest angular 
modes. Hence, rather than imaging the 21 cm bright- 
ness temperature in the target fields, reionisation sci- 
ence with the MWA will focus on statistical measures 
of the 21 cm emission. The number of independent 
measurements in a raw, multi-frequency MWA data 
cube will be ~ 10*. These samples record largely 
Galactic synchrotron and free-free emission (> 200 K) 
and extragalactic continuum sources (~ 50 K), with 
only a small contribution from the redshifted 21 cm 
signal (~ 25 mK). The primary approach planned for 
the MWA is to calculate the power spectrum of spatial 
fluctuations of the 21 cm signal. 
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Figure 1: Simulated 21 cm brightness temperature maps from simulation SI by McQuinn et al. (20071 at 
three epochs calculated under the assumption that Ts ^ T~^- The left panel is an almost fully neutral 
IGM (before reionisation at z = 12) such that the 21cm emission traces the matter density. The middle 
panel is during reionisation {xhi ~ 0.5, which occurs at 2; = 8 in this simulation): The dark regions 
represent ionized regions around galaxies. Lastly, the right panel shows the residual emission near the end 
of reionisation (z — 6). Each panel is 100 cMpc across and subtends 0.5° (a small fraction of the MWA 
field of view). 



The three-dimensional power spectrum of redshifted 
21 cm emission is the square of the Fourier trans- 
form of ST21 (Eq. [2|. It has a large spherically- 
symmetric component due to the isotropy of space, 
but has slightly more power for modes that point in 
the line-of-sight direction that results from peculiar 
velocities. The predicted components in the full three- 
dimensional power spectrum can be decomposed into 
terms with different powers of fi = k ■ n, where /x is 
the cosine of the angle between the line of sight and 
the wavevector ((Kaiser |1987| [Barkana fc Loeb.,2005 | 
McQuinn et"aL||2006[ ). Thus , 



P2I (fc) ^ P^0{k) + fl^P^2{k) + fl^P^^k). 



(3) 



During reionisation, the P^o term typically dominates, 
and the most important contrib ution to P^o is the fluc - 
tuations in the neutral fraction ( McQuinn et al.|2006 l. 

Averaging P2i{k) over angle to yield P2i(fc) loses 
some statistical information, but at the sensitivity of 
MWA this loss is minor during reionisation. With one 
year of the planned observing program, the MWA can 
constrain P2i{k) at a given redshift and across a range 
of scales sp anning 0.01 ^ k < 1 cMpc~^ with total 
S/N 14(|Morales fc Hewitt|[2004| [Bowman et al. 



2006l|McQuinn et al.|2006[|Beardsley et al.|2013p . The 

sensitivity of the MWA to the spherically-averaged 
power spectrum is shown in Figure [3] 

The properties of the power spectrum, P2i{k), dur- 
ing reionisation are dominated by the characteristics 
of the ionized bubbles (Figure [l] middle panel) that 



are determined by the ionizing radiation produced by 
the first galaxies and quasars. Figure [3] illustrates how 
A21 = (fc^P2i(fc)/27r^)^^^ contains information about 
the reionisation process. It plots A21 for the fiducial 
model of |Lidz et al.| ( |2008[ ) at several neutral frac- 
tions. The shape of the power spectrum evolves con- 
siderably as the hydrogen in the IGM becomes ion- 
ized, with power shifting from small scales initially, 
to larger scales later when large ionized bubbles have 
formed. Eventually, the amplitude decreases as hydro- 
gen reionisation completes. 

The sensitivity forecasts shown in Figure [S] utilize 
both the primary and secondary target fields in combi- 
nation, assuming that each field has been observed over 
an 8 MHz bandwidth (corresponding to Az ~ 0.5). 
The full 32 MHz instantaneous bandwidth of the MWA 
will yield four such measurements spanning a contigu- 
ous redshift range of Az ~ 2. Initial observations will 
focus on the lower redshifts (6 < 2; < 8), where the sen- 
sitivity of the instrument is highest due to the lower 
sky temperature. Observations at low redshifts will 
crucially demonstrate whether the foreground subtrac- 
tion is successful, since the quality of the foreground 
subtraction will be readily apparent at z < 6, where no 
21 cm signal should be detectable by the MWA since 



Xhi < 2% after reionisation ( |Wyithe fc Loeb|2008 1 



2.1.1 Sources of reionisation 

The 21 cm sky contains a wealth of information about 
the properties of the ionizing sources because they 
strongly affect the topology of the EoR. Reionisation 
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Figure 2: Pr ojected all-sky brightness tempera- 
ture map of de Oliveira-Costa et al.| (12008) 



150 MHz. The primary (1) and secondary (2) 
get regions for EoR observations are centred at 
a,S = 60°, -30° {£,b = 228°, -49°) and a,S ^ 
155°, -10° {£,b = 253°,-H38°), respectively The 
solid white curves indicate the 50% response power 
contour of the primary antenna tile beams for ref- 
erence and the dashed lines show the 10% power 
contour. 



scenarios dominated by stars in galaxies produce well- 
defined ionized regions, while "miniquasars" (small ac- 
creting black holes) yield more difi'use features ( |Zaroubi| 
|fc Silk 2005) and reionisation processes involving de- 
caying/annihilating dark matter particles would be very 
inhomogeneous ( [Belikov fc Hooper|[2009l ) . The 21 cm 
observations could even disti nguish between di fferent 
stellar reionisation scenarios (Lidz et al. 20081. For 



example , more massive galaxies produce larger ionized 
regions ( |Furlanetto et al.|2004| [McQuinn et al.||2007 l. 

If it is able to constrain the shape and evolution of the 
power spectrum, the MWA will reveal the nature of 
some of the first cosmological objects. 

2.1.2 The Neutral Fraction 



Existing observations have provided tantalizing hints 
about reionisation, but even more unanswered ques- 
tions. CMB observations fin d that the average red shift 
of reionisation is z = 10 ±2 ( Komatsu et al.|2009 1, and 
quasar absorption spectra imply that the Universe was 
totally ionized by z = 6, one billion years after the Big 
Bang. In addition, several studies have interpreted the 
large saturated regions at z > 6 in these absorption 
spectra as evidence for the end of reionisation (e.g.. 



Fan et al. 2006 1, but this conclusion is disputed (|Becker 
et al.|2007| |Mcsingcr 20101 



Detection of the 21 cm power spectrum will yield 
measures of xhi{z) through predictable features that 
depend on the neutral fraction. As evident in FigurejSj 
the intensity of 10 cMpc-scale fluctuations in 21 cm 
brightness temperature briefly fades when reionisation 
begins and the first galaxies ionize the overdense re- 
gions within a couple of cMpc of them. Next, the ion- 
ized regions surrounding individual galaxies overlap, 



creating ~ 10 cMpc ionized regions and an increase in 
power on these scales. Finally, all of the fiuctuations 
fade as these bubbles overlap, swallowing the neutral 
hydrogen that remains. Each of these trends can be 
used to disentangle xhi{z). 

For many possible sources that have been consid- 
ered in the literature (such as galaxies, quasars, and 
mini-halos), the amplitude of the power spectrum peaks 
at xhi « 50%. The MWA could probe the global ion- 
ization history by tracking the evolution of this ampli- 
tude and shape of the 21 cm power spectrum. Using 
these features, the volume filling factor of neutral hy- 
drogen (comparable to the mass-averaged neutral frac- 
tion) can be constrained after approximately two years 
of observations to between 0.4 < xhi < 0.75 at 2a sig- 
nificance for the redshift where the power spectrum 



achieves its maximum amplitude (Lidz et al. 2008 



While unlikely, it is possible that the IGM is highly 
ionized over the range of redshifts measured by the 
MWA. [Lidz et al.| (2008) found that a nuU detection 
with the MWA at 2; < 8 would be able to constrain 
the neutral fraction to < 4%. 

Observations of individual quasar HII regions of- 
fer an additional probe of the global neutral hydrogen 
fraction since the contrast between the ionized HII 
region and neutral IGM can provide a measurement 
of the neutral fraction at the redshift of the quasar. 
Ceil et al.| ( 2008 1 have found that the shape of the 
quasar HII regions can be recovered from foreground 
subtracted MWA maps by using the entire data cube 
to help determine the contrast between the external 
IGM and ionized interior. 



2.2 Quasar HII Regions 

The detection of large quasar-generated HII regions in 
redshifted 21 cm emission late in the epoch of reioni- 
sation may be poss ib le with the MWA ( [Wyithe et al 



2005| |Wyithe|[2008| . |Geil et al.| ( |2008[ ) have modeled 



the observation of quasar HII regions in a partially ion- 
ized IGM. Figure [4] shows a simulated example of a re- 
covered quasar HII region with the MWA after approx- 
imately one observing season. The recovery is possible 
only for quasars with positions known a priori. Al- 
though the EoR fields will not overlap the Sloan Digi- 
tal Sky Survey (SDSS) fields, where the existing known 
z ~ 6 quasars reside, the contemporaneous SkyMapper 
optical survey ( [Keller et al.||2007l ) should find around 
one quasar at z > 6.25, and around 5 quasars at z > 6 
per MWA field. 

Like all 21 cm signatures of the reionisation era, 
the detection of a quasar HII region will need to over- 
come the difficulties associated with the removal of 
bright Gala ctic and extragalactic foregrounds. |Geil| 
et al.[ {[2008) find the jrimary effect of continuum fore- 
ground removal is to reduce contrast in the image. In 
particular, contributions to the 21 cm intensity fiuc- 
tuations that have a scale length comparable to the 
frequency bandpass of the observation are removed by 
foreground subtraction. This loss of contrast, however, 
does not severely affect the ability of 21 cm observa- 
tions to measure the size and shape of the HII r egion . 
To demonstrate this quantitatively, Geil et al. ( 2008 ) 
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2.3 Non-Gaussian Statistics 
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Figure 3: MWA thermal uncertainty at 2: = 8 
with simulated 21 cm power spectra. The black 
stepped line is the la thermal uncertainty mod- 
eled for a combined observation of the primary and 
secondary EoR fields with a total of 1600 hours 
( [Beardsley eraTHMsl ). The 21 cm curves are cal- 
culated using the fiducial simulation described by 
Lidz et all (|2008|). 



The power spectrum is very 
steep at the beginning of reionisation {xhi ^ 0.8) 
and falls rapidly at the scales probed by the MWA. 
Large-scale power returns as large ionized regions 
form by xhi ~ 0.5 and then it falls again with 
decreasing neutral fraction. The MWA should be 
able to detect the power spectrum at 2: = 8 for any 
appreciable neutral fraction outside of the range 
0.95 > Xhi > 0.8. 



modeled the HII region shape as a spheroid described 
by six parameters and showed that the shape recov- 
ered following foreground removal agrees well with the 
shape derived directly from fitting using the input HII 
region model. 

Using the recovered best-fit shape of the HII re- 
gion, the global neutral fraction of hydrogen in the 
IGM can be measured directly from the contrast in in- 
tensity between regions that are within and beyond the 
HII region. However, since foreground removal low- 
ers the observed contrast between the HII region and 
the IGM, such a measurement of the neutral fraction 
requires a correction factor. The value of this correc- 
tion factor depends on the reionisation history and the 
polynomial used for foreground subtraction and can be 
accounted for using astrophysical models. 



During the reionisation epoch, fluctuations in the 21 cm 
emission are expected to be highly non-Gaussian due 
to galaxy bias and the patchy structure imparted by 
the ionized bubbles in the IGM. The power spectrum 
(or two-point correlation function), therefore, is not 
a sufficient descriptor of the properties of the field, 
and additional statistical measures can provide com- 
plementary information. Several extensions to the power 
spectrum have been proposed in this context, includ- 



ing the 21 cm-galaxy cross spectrum ( Furlanetto et al. 
Wyithe et al]|2007[ |Lidz et al.||2009[ ), direc% 



2006 



measuring the one-point probability distribution func- 
tion (PDF) of the 21 cm brightness tem perature distri- 
bution as a function of smoothing scale fMellema et al. 



2006 Ichikawa et al. 2010, Pctrovic 



Oh 2011J or its 

higher moments ( [Wyithe fc M orales 2 Q07| , the three- 
point correlation function (bispectrum), matched fil- 



ters ( Datta et al. [200^ 1, the th reshold clustering func - 
tion, and other derived metrics ( Barkana fc Loeb|2008 \ 



2.4 Foreground Subtraction 

The ability of the MWA to achieve its reionisation sci- 
ence is entirely contingent on the successful separation 
of the 21 cm signal from astrophysical foregrounds. 
Hence, developing and demonstrating the techniques 
required to subtract foregrounds is a critical research 
objective. 

2.4.1 Reference Pipeline 

The baseline foreground subtraction strategy is a multi- 
stage process consisting of the following principal stages: 

1 - Bright source subtraction. The first step in 
foreground subtraction is the removal of bright points 
sources (e.g. extragalactic AGN). This is necessary 
because the sidelobes of these sources mask the 21 cm 
signal. Recent theoretical studies have indicated that 
bright point source subtraction must be complete for 
sources with flux greater than S > 100 mjy in order 
to allow for successful subsequent Galactic continuum 
foreground subtraction (Liu et al. 2009). The bright- 
est sources in the field will be "peeled" (self-calibrate 
and subtract) from the raw visibility measurements 
( Mitchell et al.|["2008 |. Efficient techniques to remove 
point soiu-ces from grid ded and integrated maps are 
also under investigation ( Pindor et al.|[20TT Bernardi 
et al.|2011[ ). 

2 - Diffuse continuum subtraction. A variety of 
mechanisms contribute to diffuse continuum emission 
in the low-frequency radio sky including Galactic syn- 
chrotron emission. Galactic free-free emission, extra- 
galactic free-free emission, and a sea of faint extra- 
galactic point sources that are too numerous to iden- 
tify individually in maps and are, thus, said to be 
"confused" . All of these foregrounds are significantly 
stronger than the redshifted 21 cm signal and have 
considerable angular structure on the sky (see Fig- 
ure |2] and Santos et al. 2005 their Figure 5). But 
the emission mechanisms of these foregrounds produce 
smooth spectra that follow power-law profiles, making 
them possible to separate from the redshifted 21 cm 
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Figure 4: Si mulated recovery of a quasar HII region in an evolving IGM, as observed by an MWA-like 
instrument (Geil et al. 2008). Foregrounds were simulated and removed for the maps shown. The data 
cube is only a small part of a full MWA field and spans ^ 3.3 degrees on a side and 33 MHz deep. The 
three panels show different slices through the centre of the box when viewed from the front, side, and 
top. Each slice is 6 cMpc thick, which corresponds to ~ 3 MHz along the xa-axis. The extracted shape 
of the HH region is shown, based on a 9 parameter model. The quasar was assumed to have a lifetime of 
4 X 10^ years centred on 2; = 6.65 and to have produced an HH region with a radius of 34 cMpc at that 
redshift. The mass-averaged IGM neutral fraction at the redshift of the quasar is ~ 15%. 



emission by their long spectral coherence scales com- 
pared to the small, ~ 1 MHz spectral coherence of 
the reionisation signal. A variety of techniques have 
been proposed to accomplish this separation including 
polynomial fits along the spectral axis of each pixel 
in an image map or its uv-map conjugate, Minkowsky- 
functionals, and non-parametric techniques. The base- 
line approach planned for the MWA is to fit and sub- 
tract a third-order polynomial from each pixel in the 
uv-domain representations of the bright source cleaned 
maps. This technique has been demonstrated for the 
MWA in the limit of perfect instrumental calibration 



( Liu et al.|20 09. Bow man et al.|2009l |Liu et al.|2009 | 

3 - Polarized leakage subtraction. In addition to 
the foreground issues discussed thus far, gain and phase 
calibration errors and non-ideal feeds conspire to con- 
taminate all Stokes parameters further. This leakage 
leads to some portion of the complex hnear polariza- 
tion signal finding its way into the Stokes I intensity 
maps that will be used for the redshifted 21 cm mea- 
surements. By using rotation measure (RM) synthesis 
(Brentjens & de Bruyn 2005 see also § 3.31 to ex- 



ploit the Fourier relationship between the polarized 
signal and the Faraday dispersion function, it is pos- 
sible to identify components of polarized signal with 
specific Faraday depth. The three-dimensional Stokes 
I data set is first re-sampled along the frequency axis 
into the natural A^-basis for Faraday rotated emission. 
A one-dimensional Fourier transform is then applied 
along each line of sight. Any feature observed with a 
non-zero Faraday depth in the three-dimensional RM 
synthesis map will, by definition, be polarized con- 
tamination. By applying an inversion process together 
with knowledge of the instrument's specific leakage be- 
haviour, this contamination can be subtracted (Geil 
et aL|2011[ ). 

4 - Statistical template fitting and power spectrum 
estimation. Many classes of errors and residuals re- 



maining in the maps after the foreground cleaning stages 
described above can be modeled given suflicient knowl- 
edge of the instrumental response ( Morales et al.|2006 



Bowman et~aLl[2009l fPatta et al.1|2010[ [Morales et al.| 



2012 1. The statistical properties of these residuals can 



be exploited to provide an additional layer of fore- 
ground removal by fitting templates of the residual 
statistical structures during the final redshifted 21 cm 
power spectrum parameter estimation. This technique 
is now commonly employed in CMB power spectrum 
analysis to account for the angular power of faint point 
sources and the effects of gravitational lensing. 



2.4.2 Unified Parameter Estimation 

The stages of foreground subtraction outlined above 
can be performed discretely, resulting in a series of 
output products, the last of which is used as the in- 
put into power spectrum estimation and other statis- 
tical tests to extract information from the redshifted 
21 cm signal. Alternatively, they can be performed 
in a unified matrix framework that simultaneously fits 
models for the foreground contributions and the 21 
cm st atistical signal (Liu & Tegmark 2011 [Tegmark 
19971. The latter has the advantage of producing an 



optimal result that minimizes both bias (from thermal 
noise and foreground contamination) and correlations 
between parameters, and yields the smallest possible 
error bars, but at the expense of a large computational 
overhead. Advances in processing power and new ef- 
forts at sparse matrix reduction have begun to enable 
the trial application of unified parameter estimation in 
simulations and preliminary data, and it is expected 
this approach will ultimately be the preferred tech- 
nique for data analysis. 
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3 Galactic and Extragalactic 
Science 

The MWA offers several exciting science opportunities 
under the collective banner "Galactic and Extragalac- 
tic" science. The major input into this science will 
come from a deep all-sky survey over the MWA fre- 
quency range with full polarimetry and spectral reso- 
lution. This survey is a critical undertaking as it will 
also be used by the other science programs as a basis 
for an all-sky foreground model and the identification 
of calibrators. The data will also be combined with 
other new surveys, including GALFACTS and POS- 
SUM (1.4 GHz), S-PASS (2.4 GHz), and Planck (27 
GHz to 1 THz) surveys. 

The following sections describe a selection of Galac- 
tic and extragalactic experiments for which the MWA 
is well suited. These projects capitalize on three ad- 
vantages of the MWA's low observing frequencies: (1) 
non-thermal sources will be very bright; (2) Faraday 
rotation and depolarization will be very strong; and 
(3) along some sightlines free-free absorption will be 
significant. 

3.1 The Cosmic Web 



Theory suggests that the infall of baryons into the cos 
mic web res ults in high Mach number intergalacti c 
shocks (e.g 



et al.|2010 



2012 k The MWA's surface-brightness sen- 



Cen fc Ostriker 1999 Ryu et al. 20031 



When combined with seed magnetic fields, efficient 
cosmic ray acceleration is possible, resulting in steep- 
spectrum synchrotron emi ssion ([Keshet et al.| |2004| 



Pfrommer et al.|2006[[2008l|Skillman et al.|2011||Araya^ 



Melo et al. 20121. This may be a ubiquitous compo- 



nent of the extragalactic radio sky ( Keshet et al. 2004 1 



and recent simulations suggest that the majority of the 
luminous radio shocks can be found in the outskirts 



of massive galaxy clusters (Araya-Melo et al. 20121 



The MWAs centrally condensed antenna configuration 
provides the theoretical sensitivity required to make 
the maps of this relativistic cosmic web. However, the 
accurate subtraction of discrete sources requires long 
baselines and high image fidelity | |Brown|[20Tl| . This 
may limit direct to the densest portions of the web, 
though statistical detection may be possible ( [Keshet | 
et aL|2004| [Brown et aL]|2010[ ). 



3.2 Radio Relics and Clusters 

The MWA is also well suited for studying the diffuse, 
steep-spectrum synchrotron emission from clusters of 
galaxies. These objects, which include the cent rally lo- 
cated 'radio halos' and peripheral 'radio relics' ( Kemp- 
ner et al.|2004[ ), are probes of relativistic electrons and 
large-scale magnetic fields in the intra-cluster medium 
(ICM; [Ferrari et al.[[2008| [Rudnick et al.|[2009[ ). Pos- 
sible mechanisms for accelerating the cosmic ray elec- 
trons responsible for this emission are turbulent re- 



acceleration following a cluster merger ( Brunetti et al. 
2001 Petrosian|2001 1 and hadro nic collis ions bet ween 
cosmic-ray and thermal protons (Dennison 1980; Blasij 
& Colafrancesco!'1999). Low frequency observations 



of clusters can probe these different models ( jCassano 



sitivity wil enable it to detect very steep-spectrum 
halos and relics. Predictions for the number of such 
halos and relics detectable with other low frequency 
instruments are in the range of several hundred out to 
z < 0.6 ( [Cassano et al.|20 12). 

At 150 MHz, the MWA will resolve 1 Mpc struc- 
tures out to a redshift of 0.2 and 500 kpc structures out 
to 2; < 0.07. As with the detection of the cosmic web, 
confusion will be the major obstacle. Use of comple- 
mentary higher resolution observations to model and 
subtract confusing sources or will be required (Brown 
20Tl| . The TIFR GMRT Sky Survey (TGSS) that 
covers 32,000 square degrees north of declination —30° 
with angular resolution 20" presents a good compari- 
son survey to the MWA for this work. 

3.3 Faraday Tomography and Mag- 
netic Fields 

The Milky Way and most other spiral galaxies all show 
well-organized, large-scale magnetic fields, but the pro- 
cesses that amplify and maintain these fields are still 
not understood. On smaller scales, fluctuations in 
these magnetic fields regulate the dissipation of large- 
scale gas flows via a turbulent cascade. However, it 
is as yet unclear what processes generate turbulence 
in the interstellar medium (ISM), and how the subse- 
quent spectrum and outer scale of turbulence depend 
on the surrounding conditions. 

Polarized radio waves experience Faraday rotation 
when they propagate through a magnetized plasma. 
The amount of rotation depends on the square of the 
wavelength, and on the path integral of neB\\ along 
the line of sight (where We is the electron density of 
the ISM, and _B| | is the magnetic fleld strength parallel 
to the line of sight). By exploiting this phenomenon, 
the MWA will provide a new view of the large- and 
small-scale properties of otherwise unobservable mag- 
netic fields in the Milky Way. 

At the low observing frequencies of the MWA, small 
rotation measures have substantial effects on the dis- 
tribution of polarized emission from the sky. This 
emission directly traces invisible foreground magnetic 
fields, which imprint their structure on the smooth po- 
larized synchrotron emission from the Milky Way via 
foreground Faraday rotation. By imaging the polar- 
ized sky at a range of frequencies, the MWA will pro- 
vide direct measurements of the turbulent cascade in 
the local interstellar medium, and will reveal the de- 
tailed magnetic geometries of discrete large-scale struc- 
tures such as the Local Bubble, the North Polar Spur 
and the Magellanic Clouds. Using the high spectral 
resolution of the MWA, we can also apply RM synthe- 
sis techniques (see also § [2.4. 1[ ) to isolate overlapping 
RM components. An example is provided in Fig. [5| 
which shows the turbulent pattern seen at 188 MHz 
in a large region of the sky. The MWA will provide a 
unique tomographic view of magneto-ionized gas, al- 
lowing the identification of multiple emitting and ro- 
tating regions along the same line of sight. 

The MWA is expected to make the largest con- 
tributions to Galactic science by observing at mid- 
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Figure 5: Polarized intensity at Faraday depth zero at 188 MHz from the MWA development system 
(Bernard! et al. 2013). The field spans 2400 square degrees with 15 arcminute angular resolution and 



is roughly centered on the southern Galactic pole, extending to within 20 degrees of the Galactic plane. 
The radial lines are Galactic latitude and the circular arcs indicate Galactic longitude. Flux density is 
plotted on a linear color scale from to 0.2 Jy beam~^ RMSF~^ (white to black), where RMSF is the 
rotation measure spread function. Polarized maps from the MWA showing similar complex structures and 
discontinuties that often have no counterparts in total intensity will be used to determine the properties 
of both ordered and turbulent magnetic fields in the Milky Way's ISM (see Figure l6|. 



and high-latitude regions. Depolarisation is already 
severe in the inner Galactic plane at 1.4 GHz ( Reich 



et al. 1990 Wolleben et al. 2006 Haverkorn et al. 2006 1 



and it is likely that much of the inner Galactic plane 
will appear completely depolarised below 300 MHz, 
although [Bernardi et al.| (^009 ) have detected polar- 
ized emission at arcminute scales at 150 MHz. Re- 
gions at mid and high latitudes, on the other hand, 
show conside rable diffuse Galactic polarised structure 
at 350 MHz ([Haverkorn et al.||2003b[ |de Bruyn et al 



2006 



Schnitzeler et al. [2009 1. The high sensitivity 



of the MWA to extended structure will enable it to 
map the diffuse polarised emission down to very low 
intensities, thus making visible weak Galactic mag- 
netic fields in the local interstellar medium and Galac- 
tic extended disk and halo that are undetectable at 
higher frequencies. RM synthesis across the entire 
MWA frequency band will result in a RM resolution 
of 0.3 rad m~^, which gives unprecedented detail in 
RM structure, allowing detailed studies of e.g. inter- 
stellar magnetised turbulence, or magnetism in old 
supernova remnants. Increasing depolarisation as a 
function of decreasing frequency will give independent 
measurements of the amount and scale of fluctuations 



in the magneto-ionised medium (Sokoloff et al. 1998 
Gaensler et al.|2001| . 



A powerful new additional probe of the magne- 
tised ISM comes from the polarisation gradient, in 
which the spatial derivatives of Stokes Q and U are 
combined to reveal locations on the sky at which the 
electron dens ity or magnetic field strength show sharp 
fluctuations (Gaensler et al. 2011). Application of the 



polarisation gradient at higher frequencies reveals a 
complex tangled network of narrow filaments (see Fig- 
ure [gJ . These structures are a superset of the "polari- 
sation canals" seen in many radio polarisation images 



dHaverkorn fc Heitsch|2004l|Fletcher fc Shukurov|2007 
Schnitzeler et al.||2009[ ),"and directly reveal the loca- 
tions of shocks, cusps and shear in the turbulent ISM 
Calculation of the polarisation gradient for magne- 



tohydrodynamic simulations of interstellar turbulence 
shows that the intensity, alignment and morphology 
of these filaments depend sensitively on input parame- 
ters such as the sonic Mach number and magnetic field 
strength ( Burkhart et al.||2012 l. By applying statisti- 
cal metrics such as the genus and kurtosis, one can 
directly compare observations of the polarisation gra- 
dient with a grid of simulations, and thereby extract 
quantitative information on the otherwise unobserv- 
able turbulent parameters of the ISM. 

Two effects distinguish the applicability of this tech- 
nique for MWA data. First, the "polarisation horizon" 
( [Uyaniker et al.|2003| ) at the MWA's low observing fre- 
quencies is very nearby (~ 100 pc, depending on ob- 
serving frequency and Galactic coordinates; |Bernard"i| 
et al. |2013[ ), meaning that we can probe small-scale 
local features. Second, the change in electron den- 
sity or magnetic field strength needed to produce a 
fluctuation in Q or f/ is very small (since the amount 
of Faraday rotation is proportional to A^). For these 
two reasons, application of the polarisation gradient 
to MWA data should provide a detailed view of tur- 
bulence in the solar neighbourhood, and will be able 
to extend studies of the global turbulent cascade to 
smaller amplitudes and smaller angular scale than has 
been possible from previous observations. 

Rotation measures from extragalactic polarised point 
sources can be used to probe the intervening medium 
between the source and the observer, where gener- 
ally the Galactic RM component dominates. Little 
is known about the polarisation percentages of extra- 
galactic point sources at low frequencies. At 1.4 GHz, 
the expected number of linearly polarised sources down 
to a 3(7 s ensitivity of 1 mjy is about 1.2 x 10** per 
steradian ( Beck fc Gaensler|2004 1 . The number of po- 
larised sources may decrease at low frequencies due 
to increased inter nal depolarisation. At 1.4 GHz, 
Mesa et al. ( 2002 1 found a median polarisation frac- 
tion of 2.2% for sour ces above 80 mJy in the NVSS 
and Taylor et al. ( 2009| ) have compiled a catalog from 
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Figure 6: Gradient image of linear polarization, |VP|, for an 18 deg^ region of the Southern Galactic 
Plane Survey ( [Gaensler et al.|2011 ) that has been used to study magnetized turbulence in the interstellar 
medium. Similar statistical techniques will be applied to MWA polarized maps to investigate magnetic 
field strength, Mach number, Reynolds number, and othe properties in the interstellar medium. 



the NVSS of over 30,000 polarised sources detected 
at better than 80" and selected for polarisation frac- 
tion larger than 0.5%. The average polarisation frac- 
tion of sources in this catalog is ~ 7%. At 350 MIfz 
[Schnitzeler et ah] ( |2009[ ) found an average polarisa- 
tion fraction of 3 % for polarized sources b rig hter than 
3 mjy (see also Haverkorn et al. 2003c I. Bernardi 
( [20T3I ), have surveyed 2400 square degrees at 



et al. 



188 MHz using the 32-element MWA development sys- 
tem. They achieved noise levels of 15 mJy in polarisa- 
tion and found only one polarized source, PMN J0351- 
2744, which had a polarization fraction of 1.2% and 
an RM of -1-34 rad m"^ The source was detected at 
20 a. They also placed an upper limit of ~ 1% on the 
polarisation fraction of over one hundred other bright 
sources above 4 Jy and detected diffuse polarised emis- 
sion from the Galaxy with a peak level of 13 K. 

Beam depolarization is likely to be partially re- 
sponsible for the lack of detected sources by |Bernardi| 
et al. (20131 since the MWA demonstrator had only 



15 arcmin resolution. Using 1% polarization fraction 
as a lower limi t, along with sou rce flux counts from 
the 6C survey ( Hales et al.||1988 l and the ratio of po- 
lar ised sources to unpol arized sources of 14% found 
by Banfield et al. (20111, we estimate the number of 



sources that will be detectable by the MWA assuming 
a fiducial sensitivity of 1 mJy to be approximately ~ 1 
source per square degree. These extragalactic source 
RMs will provid e a "grid" ( |Beck fc Gaensl er 2004; Tay-| 
lor et al.|[2009t that can be used to model large-scale 
magnetic fields in the Milky Way and, if deeper obser- 
vations can provide a sufficiently dense grid, in nearby 
galaxies. Provided that polarisation calibration is ac- 
curate enough, the MWA will observe point source RM 
values with superb precision. 

3.4 The Magellanic Clouds 

The Magellanic System, consisting of the Large Mag- 
ellanic Cloud (LMC), Small Magellanic Cloud (SMC), 



and other components, such as the Magellanic Stream 
and the Bridge and Leading Arm, is the Milky Way's 
most significant neighbour. The Clouds and the Milky 
Way are engaged in a three-body interaction allow- 
ing us to study at close hand the effect of dynamical 
interaction processes on the physical structure, star- 
formation, magnetic fields and chemistry of such galax- 
ies. The location of the Clouds in the deep southern 
sky makes them ideal candidates to observe with the 
MWA in preparation for much more detailed studies 
at similar frequencies with the SKA. 

A comparison of radio and far-infared emission in 
the Magellanic Clouds and other nearby galaxies al- 
lows for a more detailed study of the correlation of 
cosmic rays, magnetic fields and star-formation than 
is possible for distant galaxies. At GHz frequencies, 
ATCA and IRAS studies allowed Hughes et al. (2006) 
to find that the radio/far-infrared (FIR) relation breaks 
down at ~ 50 pc (3 arcmin). Residual correlation at 
smaller scales remains due to the thermal component 
of the radio emission. Low-frequency observations of 
6C galaxies have shown that the radio/FIR correlation 
remains good at 151 MHz (Fitt et al. 1988). Such ob- 
servations are useful for further separating thermal and 
non-thermal sources of emission, for the study of syn- 
chrotron spectral aging and cosmic ray diffusion, and 
for allowing k-corrections in cosmological studies. 

The Magellanic Clouds are among only a handful 
of galaxies where the magnetic field can be studied 
via both Faraday rotation of the polarized emission 
from background sources (sensitive to the line of sight 
field) and the polarized component of the diffuse emis- 
sion (sensitive to the field direction in the plane of 
the sky). Studies of the Clouds at GHz frequencies 



dHaynes et al.||1986[ [Gaensler et al.||2005a| |Mao et~al 



2008 1 have been useful in determining field strengths 



and morphology. Low-frequency observations will con- 
tribute by locating regions, perhaps associated with 
wind-blown superbubbles, of warm ionized gas unseen 
in X-rays, and by allowing more detailed studies of 
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Faraday structures. The MWA will investigate the ex- 
istence of a possible pan-Magellanic magnetic field via 
the old relativistic electron population which will have 
travelled far from its source. 



the full energy budget of the Galaxy, including the role 
of supernovae in the production of energetic cosmic 
rays, as well as the diffusion and aging of the cosmic 
rays. 



3.5 Cosmic Ray Mapping 

One of the most exciting possibilities of the MWA sur- 
vey at the lowest frequencies is the mapping of the 
three-dimensional cosmic ray emissivity of the Milky 
Way. Within the MWA frequency range below ap- 
proximately 200 MHz, Galactic HII regions turn op- 
tically thick, blocking emission from the synchrotron 
background of the Galaxy. Typical optical depths of 
HII regions are expected to be r « 2 to 5. These 
HII regions, therefore, appear in absorption relative to 
the surroundings and the residual surface brightness 
at these positions allows a direct measurement of the 
integrated emissivity of cosmic-ray electrons out to the 
HII region. 



Nord et al. (20061 used 74 MHz VLA observations 



of HII regions to detect such absorption against the 
Galactic synchrotron background. They found 92 re- 
gions with absorption features, 42 of which have mea- 
sured distances. While such features have been seen 
since the 1950s, previous observations had poor reso- 
lution. The 10 arcmin VLA beam allowed for better 
spatial resolution while still remaining sensitive to low 
surface brightness. A number of surveys at higher fre- 
quencies have been conducted. Paladini et al. (20031 



3.6 Galactic Supernova Remnants 

Statistical studies of supernova rates suggest that there 
should be ~ 1000 — 2000 supernova remnants (SNRs) 
in our Galaxy ( |Li et al.|[l991| [Tammann et al.||1994[ ), 
in stark contrast to the ^ 300 S NRs currently known 
( |Green|2004"| [Brogan et al.|2006[ ). This deficit is likely 
due to the observational selection effects, which dis- 
criminate against the identification both of old, faint, 
large SNRs, and young, small SNRs. Discovery of 
these 'missing' SNRs, and hence a characterisation of 
the full SNR population, is crucial for understanding 
the production and energy density of Galactic cosmic 
rays and the overall energy budget of the interstellar 
medium. 

Recent studies have demonstrated that low-frequency 
interferometric maps are a very efficient way of iden- 
tifying new SNRs, because of the comparatively steep 
radio spectra of t hese system s. Mos t notable has been 
the survey of Brogan et al.| ( [20061 ), shown in Figj?] 



have compiled a catalog of 1442 diffuse and compact 
HII regions based on 24 prior surveys, although a num- 
ber of regions are omitted in their catalog due to their 
emphasis on regions of order 10 arcmin (for an ex- 
ample of more recent identifications, see [Bania et al.| 



20101. Approximately 800 of these regions have radio 
line velocity measurements providing distances. |Pala-| 
dini et al. ( 2004 1 report on the distribution of the HII 



which tripled the number of known SNRs in a 40 deg 
patch of the Galactic plane. The MWA will have a 
unique view of virtually all the inner Galaxy. The 
high-quality images that we expect to derive, across 
a range of frequencies, will provide a superb data set 
from which to identify hundreds of previously uniden- 
tified SNRs. Additionally, attempts to uncover old, 
faint sources will be greatly enhanced by development 
of new detection algorithms which we plan to imple- 
ment on MWA data (e.g., circle Hough transforms, see 
[Hollitt fc Johnston-Hollitt|2012[ ). These sources wiU be 
of interest both as a statistical ensemble, and because 
among them will be individual objects with unique or 
extreme properties. 



3.7 Radio Recombination Lines 



regions, most of which have Galactic latitude within 
±2 degrees of the plane and are distributed throughout 
the disk, with the spiral arm structure of the Galaxy 
evident in their distribution. 

Because many Galactic HII regions have well-constrainedLow-frequency radio recombination lines (RRLs) are 
distance estimates, MWA measurements of the syn- a sensitive tool for probing interstellar plasma. Be- 
chrotron absorption features will allow for a three- low 100 MHz, carbon RRLs have been detected, but 
dimensional ma p of the relativis tic gas content in the only in absorption, indicating an association with cold 

(< 50 K) neutral gas or molecular clouds. At fre- 
quencies above 100 MHz, both carbon and hydrogen 
recombination lines have been observed, with the hy- 
drogen lines likely associated with normal HII regions 
as at higher frequencies. Both carbon and hydrogen 
lines have been found in observations along the Galac- 
tic plane as far as ±4 degrees off the plane in Galactic 
latitude. 

Measurements of RRLs in the MWA frequency range 
will provide an opportunity to study both emission 
and absorption processes. They will enable new stud- 
ies of diffuse photodissociation zones that lie between 
molecular clouds and the general ISM and yield esti- 
mates for electron temperature and density in these 
regions. Carbon RRLs at 7 5 MHz exhibit a line width 
range from 5 to 50 km s~^ (Erickson et al. [1995), cor- 
responding to 1-12 kHz. Although, the MWA spectral 



Milky Way (e.g. |Sun et al.|2008[ ). For HII regions with 
kinematic distance degeneracies, MWA detections may 
be able to break the ambiguities. Since the MWA will 
have a view of the inner and outer Galaxy, it can pro- 
vide the first accounting of the overall distribution of 
Galactic cosmic rays and of the large-scale magnetic 
fields that cause them to emit. In addition, due to 
its ability to measure the absorption in HII regions 
as function of frequency, the MWA should be able to 
yield the cosmic ray electron emissivity spectrum in 
the directio n of the HII regions. As noted by |Nord| 
et al. ( 2006 1, multi-frequency observations should lead. 



therefore, to the mapping of the emissivity spectrum 
of cosmic rays. 

Combined with the MWA's expected abilitity to 
detect and characterise supernova remnants (described 
below) these measurements will open new probes into 




Figure 7: A 330-MHz panorama of part of the Galactic first quadrant, derived by smoothing VLA data 
from Brogan et al. (20061 to 3 arcmin resolution. At these low observing frequencies, most of the bright 
sources correspond to non-thermal emission from SNRs. The MWA will be able to carry out such obser- 
vations over almost the entire inner Galaxy. 



resolution is not ideal for such narrow lines, single- 
pixel RRL observations have previously been used with 
much success ( Alves et al.l 2012 1. Sensitivity can also 



be regained by stacking different transitions within the 
observing band. 



4 Time-domain Astrophysics 

Systematic exploration of the time domain is becom- 
ing a dominant frontier in astronomy, and significant 
resources have been deployed across the electromag- 
netic spectrum in order to probe and understand tran- 
sient and time-variable phenomena. Radio transients 
vary on time scales ranging from less than one sec- 
ond to many days. The sources are compact objects 
or energetic/ explosive events, and each type is being 
investigated at the forefront of astrophysics. 

Radio observations can provide keys to understand- 
ing the explosive and transient Universe by directly 
measuring the effects of magnetic fields and non-thermal 



processes that drive the outbursts (e.g., Cordes et al. 
[2004] ). Even more intriguingly, synergies with multi- 
messenger astronomy (correlating electromagnetic tran- 
sients with events detected in gravitational waves or 
neutrinos) have the potential to unlock substantial new 
physics. The MWA spans the upper frequencies asso- 
ciated with coherent radio sources and the lower fre- 
quencies of non-thermal magnetohydrodynamic pro- 
cesses. This makes observations of transients with the 
MWA particularly useful for understanding the non- 
equilibrium processes that drive dynamic astrophysical 
systems. 

The potential sources of radio transient emission 
fall into five broad categories: coronal emission from 
nearby stars/substellar objects, emission from com- 
pact objects such as neutron stars and accreting black 
holes, explosive events such as gamma-ray bursts (GRBs) 
and radio supernovae, planetary emission from nearby 
systems, and new phenomena discovered through the 
opening of new regions of observational phase space. 



Of these, the first three are known sources of transient 
radio emission (at least at some frequencies), while the 
last two categories constitute more speculative targets. 
However, there are firm theoretical expectations that 
these latter objects should give rise to transient ra- 
dio emission, and our searches will be far deeper than 
previous surveys. 

First, we discuss in detail some of the expected 
physical constraints to be gained from the MWA ob- 
servations of known sources of transients (§ |4.1[ ). We 
then explore some of the more speculative science tar- 
gets that we hope to explore with the new phase space 
opened up by the MWA 



4.21 



We describe the new 
capabilities afforded by the MWA in § |4.3| and in § |4.4| 
we discuss observations of pulsars, which can be both 
steady and transient sources. 

4.1 Known Transient Targets 

4.1.1 Low-Mass Stars and Brown Dwarfs 

Magnetized plasmas in stellar coronae (the outer lay- 
ers of a star's magenetosphere) produce strong radio 
emission through non-thermal, collective, and/or co- 
herent processes. This emission refiects local proper- 
ties such as density and magnetic field strength, and 
can be used to understand magnetic field generation, 
magnetospheric structure, angular momentum evolu- 
tion, magnetic flaring mechanisms and accretion pro- 
cesses. Furthermore, magnetic fields provide one of the 
few indirect probes of stellar interiors. 

While radio studies of late-t ype stars an d active 
binaries are well-established (see 'Giidol "2002 for a re- 
view) , radio observations of the lowest mass stars and 
brown dwarfs — objects incap able of sustained core 
hydrogen fusion (Kumar 1962[ ) — have only recently 
been demonstrated, and comprise a relatively new tool 
in this active field of research. M dwarfs with masses 
0.1-0.5 M© are a dominant stellar component of the 
Galaxy, comprising ~70% of all stars. They are also 
known to exhibit high levels of magnetic activity, with 
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Figure 8: Phase space for radio transients observ- 
able with the MWA at a nominal frequency of 
200 M Hz, updated and adapted from |Cordes et al.] 
( 2004[ ). We plot the product of the observed peak 
flux density S'poak and the square of the distance 
(like a luminosity) against the product of the 
emission frequency v and the transient duration 
W . In the Rayleigh- Jeans approximation, these 
quantities are directly proportional and related to 
the brightness temperature T (given by the di- 
agonal lines); we use a brightness temperature of 
10^^ K (thick diagonal line) to divide coherent and 
incoherent processes (Readhead 1994). The red 



lines show the predicted sensitivity of the MWA, 
assuming a 50mJy source can be detected in 8 s. 
The vertical red lines give the timescales for indi- 
vidual snapshot observations (8 s), and short- and 
long-term surveys (5 m and 1 week). The red diag- 
onal lines give the S^cakD^ limits corresponding to 
distances of 10 pc (appropriate for sources such as 
low-mass stars, brown dwarfs, and planets) , 1 kpc 
(local Galactic sources), 10 kpc (the Galactic Cen- 
ter), IMpc (local group) and 1 Gpc (extragalactic 
sources); the dotted lines show an estimated confu- 
sion limit of 10 mJy, although we can search below 
this limit using image subtraction or for polarized 
emission 



strong surface magnetic fields (~2-4kG, compared to 
^10 G on the Sun;|Saar & Linsk^ WS ; Johns-KruU &| 
|Valenti|[T996l ) and substantial filling factors (50-80% 
of the surfaces of M dwarfs are cov ered by magnetic 
fields; Johns-KruU & Valenti [1996 1. Magnetic radio 



emission from M dwarfs has been observed in many 
nearby systems, both fiaring and quiescent, predomi- 
nantly in the 1-10 GHz band. These studies have con- 
firmed the existence of large-scale, organized mag netic 
fields around stars with fully convective interiors ( Do- 
|nati et ar||2006[ ), despite t he loss of the aQ, dynamo 
that powers the Solar field ( |Parker|1955l ). 

For even lower mass stars and brown dwarfs, the 
detection of both quiescent and flaring radio emission 
was initially a surprise (Bcrgcr ct al. 2001), as the flux 
from these sources is orders of magnitude greater than 
that expected from the Giidel-Benz relatiorj^ Several 
late-type sources have now been detected at GHz fre- 
quencies (Fig. [8|, in many cases despite t he absence 
of optical or X-ray emission (see Berger 2006 for a 
review). Indeed, while Ha and X-ray luminosities rel- 
ative to bolometric luminosities appear to decline with 
spectral type, relative radio emission inexorably in- 
creases up to the substellar regime (see Burgasser fc 



Putman|[2005l |Berger||2006"l [Audard et ai.||2007p . Un 



derstanding why radio emission is so prominent among 
these very low mass objects, and whether such emis- 
sion correlates with age, rotation or other physical pa- 
rameters, remain outstanding questions. 

Radio flares in particular demonstrate that this 
emission is associated with magnetic events, occurring 
over timescales of microseconds to minutes (e.g., Os- 



ten & Bastian 20061. The brightness temperatures 



of bursting emissions frequently exceed 10 K, indi- 
cating coherent processes as conflrmed by high cir- 
cular polarizations (30-100%). Yet radio variability 
from these sources is remarkably diverse, from strong 
(30 mJy peak flux) , narrowband, fully polarized, burst- 
ing emission with frequency drifts (Burgasser & Put- 
'man'2005 |); to perio dic bursts suggestive of pulsar-like 
beaming ( [Hallinanet al. 2007^ and Fig . [8|; to order-of- 
magn itude variations spanning years (Antonova et al. 
2007). It has even been suggested that the observed 
"quiescent" emission from these sources may be the 
result of low-le vel, sustained electron cyclotron maser 
emission (e.g., [Hallinan et al. 



2007). The physical 



origin of transient radio bursts from low-mass stars, 
their incidence rates (duty cycle estimates range from 
1-30%) and their spectral characteristics are either 
poorly constrained or entirely unknown. These ques- 
tions have a direct impact on whether low mass stars 
comprise a signiflcant fraction of radio transients in 
gener al (e.g., |Kulkarru fc Phinney|2005 Kaplan et al. 



2008) 



The MWA has the potential to contribute substan- 
tially to our understanding of both quiescent and flar- 
ing magnetic emission from low mass stars and brown 
dwarfs. By sampling the 80-300 MHz band, MWA ob- 
servations can be used to constrain the peak frequency 
(!/pk) and spectral indices about !/pk, enabling determi- 



correlation between radio fiux and X-ray luminosity 
that extends over several orders of magnitude; see |Giidel fc| 
|Benz|1993| 
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Figure 9: P-P diagram showing the distribution of 
the known pulsar population. The diamonds show 
the magnetars (SGRs and AXPs) and the stars 
the RRATs; binary systems are circled. The di- 
agonal lines show loci of constant dipole magnetic 
field (solid) and spin-down age (dotted) , while the 
dot-dashed line is an approximate "death-line" for 
pulsar activity. The MWA will be uniquely sensi- 
tive to the slower, transient population while main- 
taining sensitivity to a considerable fraction of the 
millisecond pulsar population in the lower left. 

nation of coronal field strengths (currently estimated 
at 10-100 G; Burgass er fc Putm an 2005: Bcrger 2006|, 
electron densities and the fun damental nature of the 
emission. Audard et al. (20071 argue that radio emis- 



sion from very low mass stars may in fact peak in the 
~100MHz range if it arises from electron cyclotron 
emission from weak field s (z/pk = Vc = 30— 300 MHz for 
10-100 G fields; see also |Osten et al.|2006| ). The sensi- 
tivity of the MWA is essential for this work, given the 
intrinsic faintness, but variable nature, of radio emis- 
sion from low-mass stars; estimates range over 10 ^Jy 
to 100 mjy at 100 MHz at a distance of 10 pc, depend- 
ing on the emission mechanism. Similarly, the MWA's 
wide-field capabilities are necessary to build up statis- 
tically robust samples, particularly given the roughly 
uniform distribution of nearby low-mass stars in the 
vicinity of the Sun. High spectral resolution observa- 
tions enabled by the MWA will permit studies of field 
dynamics and remote determinations of physical pa- 
rameters such as the heights of magnetic reconnection 
regions, electron speeds, and the density profile of the 
stellar coronae. 



4.1.2 Magnetars 

While a variety of neutron stars can produce transient 
or steady radio emission (discussed in more detail in 



4.4 1, some subclasses of neutron stars can give rise to 
extremely energetic, explosive phe nomena. With mag- 
netic fields of ~10^* G, magnetars (Duncan & Thomp- 



son||1992 Thompson fc Duncan|[l995 1 undergo catas- 
trophic reconnection events known as "giant fiares" 
that release > 10*'' erg. With each flare the magne- 
tars release blobs of relativistic plasma (Frail et al. 



1999||Gaensler et al.|2005b| [Cameron et al.|2005| ) with 
peak fluxes of hundreds of mJy at frequencies from 
200 MHz up to 100 GHz (Fig. |8|. Such events are 
detected through high-energy satellites, but the radio 
monitoring provides vital diagnostics of the flare prop- 
erties and energetics. 

Beyond this, some magnetars unexpectedly have 
steady radio emission that is modulated at the X-ray 
pulse period (2-12 s). Prior searches for steady-state 
radio e mission from magnet ars had not been success- 
ful (e.g., [Kaplan et al.|2002"l ), but [Halpern et al.[p005 l 
found a radio source coincident with the magnetar 
XTE J1810-197 (also [Camilo et aL][2006[ and oth- 
ers). Monitoring the magnetars through both indi- 
vidual pulses and phase-averaged emission, looking for 
both giant flares and the onset of bright phases, will 
help understand the duty cycles of the radio-bright 
phases and hence unravel the puzzles that they present. 



4.1.3 X-ray Binaries 

Relativistic outflows or jets are a widespread yet tran- 
sient consequence of accretion onto compact objects. 
Radio jets are seen for black-hole binaries (BHBs) and 
some neutron-star binaries (NSBs). These sources are 
collectively called "microquasars" (Fig. |8|, since they 
exhibit scaled-down and rapidly evolving versions of 
the jets in active galactic nuclei. Microquasars ex- 
hibit several types of jets with different ranges of radio 
power and ejection timescale. Their study is crucially 
dependent on the correlated variations in radio and 
X-ray properties. 

The most common type of jet in microquasars is 
the quasi-steady jets associated with the X-ray hard 
state ( [Fender et"aL][2004[ [McClintock fc Remillard[ 
[2006') that can persist for months in some sources. 
Steady jets are estimated to be mildly relativistic (i.e. 
~0.7c; see Fender'20061), but there is uncertainty about 
the baryonic content and hence the scale of the mass 
outflow and the fraction of accretion energy needed to 
power such jets. 

Impulsive ejections are also produced by BHBs and 
some NSBs and the brightest of these are seen as bipo- 
lar radio jets that move with velocities as high as 0.98c 
( [Mirabel fc Rodriguez 1999 1. Imp ulsive jets are associ - 
ated with X-ray state transitions ( Fender et al.|2004 |. 
In two BHBs superluminal radio ejections, coinciding 
with major X-ray fiares, led to the appearance of X-ray 
jets (more than one year later), that were imaged with 
Chandra as the jets impacted the local ISM (Tomsick| 
et al.|2003| [Corbel et al._2005) . Smaller and more fre- 



quent ejections in BHBs are seen as correlated X-ray 
and radio fiares with ~10min timescales and radio- 



band delays ( Eikenberry et al. 2000 Wilms et al. 2007 1 



The available X-ray data usually reveal recurrent and 
abrupt changes in fiux and spectral components that 
associate the small jets with instability cycles in the 
inner accretion disk. 

Study of the connection between disks and jets is 
fundamentally based on correlated X-ray and radio be- 
haviour, as outlined above. During the last decade, 
the primary limitation for this enterprise has been the 
absence of dedicated radio instruments that can part- 
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ner with X-ray facilities, but the MWA wiU radically 
change this landscape. Additional partners for jet in- 
vestigations arise with Fermi (10 keV to 300 GeV), 
and the recent expansions of TeV Observatories (e.g., 
HESS and VERITAS; 50 GeV to 50 TeV). 

Recent measurements with the GMRT (150-1450 
MHz) in India have shown variable radio emission in 



GRS 1915+10 5 dlshwara-Chandra et al.|2005[ see Fig.[8|l 
and Sco X-1 ( [Pandey et al.||2005p , with spectra that 
do not exhibit self-absorption cutoff at low frequencies. 
This is interpreted as emission from the outer regions 
of a compact jet and from an optically-thin impulsive 
jet, respectively (also see Pandey et al.|[2006 20071. 

Overall, the MWA has the capability to monitor 
many known sources for which X-ray activity is de- 
tected, establishing outburst statistics and extending 
the spectra down to low frequencies where the total 
energetics are constrained. These observations can 
help to measure the burst luminosity function and the 
X-ray/radio correlation, and investigate their depen- 
dences on various intrinsic source parameters such as 
mass and X-ray luminosity. The MWA will provide 
data archives for the two microquasars that are reli- 
able sources of jets: SS 433 and GRS 1915-1-105, as well 
as the BHBs and NSBs with the fastest outburst re- 
currence time (1-2 years: 4U 1630-47, and GX339-4, 
Aql X-1, 4U 1608-52, and the rapid burster). These 
monitoring observations will also find new outbursts 
of known sources, acting as a trigger for observations 
across the electromagnetic spectrum. Finally, the MWA 
has the potential to discover new sources: the large 
majority of the known BHBs and candidates are tran- 
sient sources that were discovered because of their X- 
ray outbursts | |Remillard fc McCUntock|2006j ). 



4.2 Serendipitous Transient Sources 

Beyond the source categories discussed above, there 
are a number of object classes which we might hope 
to see (see [Cordes et al.|2004| and [Fender et al.|2006| ) . 
These vary from the local to the cosmological, and the 
population predictions also span a wide range. Cur- 
rently, low-frequency transient radio emission has been 
proposed but not observed for these objects, and there 
are only limited observational constraints. Such sources 
are among the most exciting possibilities for the MWA 
transient surveys (§ 4.3 1. Those surveys cover a large 
area of unexplored phase space — orders of magni- 
tude in several different parameters — and provide a 
real opportunity to discover unique sources. 



4.2.1 Extra-Solar Planets 

In the Solar System, Jupiter is a strong source of low- 
frequency radio emission produced by cyclotron maser 
processes in its magnetosphere ( |Gurnett]|1974[ |Wu fc| 



Lee 



of t 



19791. Observations have established that many 



le Jovian planets discovered around other stars 



(e.g., Marcy et al. 


2005 1 have magnetic fields ( Shkol- 


nik et al. 


20031, a 


though the field strengths are not 



known. By analogy, many of those planets should have 
similar emission. This would then be a novel way of di- 
rectly detecting the planetary companions, something 



that h as only been possible in a few special cases oth- 
erwise ( Gharbonneau et al.|2002 1 . This would allow an 
unprecedented probe of the planetary systems, such as 
estimating properties of the planetary magnetosphere, 
measuring its rotation rate, and observing the interac 
tion between the star and planet ( [Lazio et al. 20041 



Griefimeier et al. (20071 discuss theoretical expecta- 
tions for radio emission from known extrasolar planets. 
Most of the sources are predicted to emit at the micro- 
Jansky level between 10-200 MHz. A handful of plan- 
ets are expected to emit at the mjy level in the MWA 
bandpass. Their results are heavily model-dependent 
and MWA observations of the currently known exo- 
planet population will help constrain further models. 
There have been a number of searches for emis- 



sion from such planets (e.g.. 


Farrell et al. 1999 


Bas- 


tian et al.||2000| |Zarka et al. 


|2001| ILazio et al.||2004| 


Lazio & Farrell|2007 1. There are no detections to date. 



although the results are not very constraining. The 
emission is predicted to be very sporadic and bursty, 
and it is difficult to search effectively for such emission 
with pointed instruments like the VLA. Instead, new 
instruments like the MWA are li kely to be the be st 
way to detect extra-solar planets ( Lazio et al.||2004 K 



4.2.2 Explosive Events 

Gamma Ray Bursts and radio supernovae may both 
produce short and long duration transient radio sig- 
nals. The afterglow emission from each type will pro- 
vide signals at the upper frequencies accessible to the 
MWA that are delayed from the initial explosion (due 
both to plasma in the ISM and to changing opaci- 
ties) and slowly brighten over a few weeks to months 
(Morales et al. 20051. These afterglows can be im- 



portant diagnostics for explosion energetics (Soderberg| 
1^81 2007, 2010) and can help understand particle ac- 
celeration in the nearby Universe ( [Chakraborti et al. 
|2011 ). While the long timescales at low frequencies 
make the MWA unfavorable for discovery of such events, 
the rarity of the more extreme ones mean that the very 
large FOV may compensate 

GRBs, supernovae, and other explosions may also 
produce prompt pulses of coherent emission during the 
initial explosio n ( [Colgate fc Noerdlinger|[T971[ [Nakar[ 
fc Piran|2011 1. The theorized prompt signals are pro- 



duced by coherent emission near the explosive shock. 



Proposed mechanisms include current oscillations ( Uso 



fc Katz 2000), synchrotron maser activity behind the 



shock ( [Sagiv fc Wa xman 2 002| ) , or collective plasma 
mode s as seen in solar bursts ([Moortgat fc Kuijpers 



2005[). Models by [Hansen fc Lyutikov| ( [2001| ) and |Benz 



fc Paesoid[ ( [1998| l also predict detectable emission, de- 



pending on the plasma environment and dynamics of 
the GRB. Some of these signals will be accompanied 
by bursts of radiation at other wavelengths, while some 
may offer the potential for multi-messenger astron- 
omy by correlating radio and gravitational-wave sig- 
nals ( [Nakar fc Piran|2011| [Palenzuela et al.[[2010[ ). 

Detection of such prompt emission from GRBs of- 
fers a number of unique insights into the properties of 
the burst s, the inter galactic medium, and cosmology 
( |Morales et al.|2005 1. First, the details of the emission 
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that we see depend heavily on the specific model for 
the shock structure, as discussed above. Second, we 
can use the emission as a probe of the media through 
which it propagates: we will see delays and dispersion 
from the local plasma, the host galaxy's ISM, the in- 
tergalactic medium (details of which can depend on 



cosmic reionisation), and the Milky Way's ISM (Inoue 
20041). 



Past searches for prompt radio emission from GRBs 
( Dessenne et al.|[l996 Benz fc Paesold||1998 Balsano 



1999 1 were all unsuccessful. Other projects attempted 
to perform blind surveys for radio transient s ( |Amy| 
et aI.|1989[|Baird et al.|1976[|Katz et al.|2003| . All of 

these attempts have suffered from constraints inherent 
in narrowband receivers and simple transient identi- 
fication sy stems. More modern projects like CASE 
and ETA ( [Morales et"aL][2005} [Patterson et al.||2008 l 
were designed to correct many of these problems, with 
digital baseband recording and interference rejection 
through de-dispersion, but the collecting areas are still 
very small. 

4.2.3 New Opportunities 

Above we discussed two classes of sources that are 
expected to emit low-frequency radio transients, and 
where detection of such transients would enable us to 
make significant advances in our understanding of such 
objects, but where we do not know the source prop- 
erties enough to consider them "guaranteed" science. 
Beyond even those classes are additional objects that 
may be found with blind searches (e.g., neutron star 
binary inspiral events that may produce radio precur- 
sors; [Hansen fc Lyutikov[[2001[ ) . 

To date, systematic transient searches with radio 
telescopes have yielded some results, but they have 



largely been difficult to interpret (e.g., Matsumura et al. 
[2007[p009| , due to poor statistics, limited localization, 
and a long time before follow-up. However, there have 
been successes that, while enigmatic, hint at the po- 
tential offered by the MWA. We do not give a full re- 
view here, but point to particular detections that hold 
promise for the MWA. 

While searching for transient events at the Galactic 
Centre (GC), Hyman et al. (20021 found two sources. 



one ve ry faint (sub-mjy) transient source. Jaeger et al. 
('2012 ) did repeated observations at a frequency of 325 MHz 
(closer to the MWA band) and did identify one signif- 
icant transient of unknown origin, although they con- 
sider several possibilities including a distant flare star. 
Only recently have surveys moved beyond examina- 
tions of archival data (e.g.. Bannister et al.|2011 ), with 
a wide range of flux limits and timescales, to system- 



atic examinations of th e variable sky (e.g., Lenc et al. 
Croft et al.[[20Tl| ). It is becoming clear that ra- 



2008 



dio transients, are rare but promissing targets that will 
lead to advances in a number of areas of astrophysics. 



4.3 Planned Searches 

We will conduct a continuous search for transient ra- 
dio emission in each image produced by the MWA's 
real-time system, looking for changes from the previ- 
ously observed sky on timescales ranging from 8 s to 
months (eventually using the high-quality integrated 
data from the sky surveys). This search will be more 
than six orders of magnitude more sensitive than previ- 
ous transient surveys in the band and will cover a much 
broader range of frequencies, transient durations, and 
dispersion delays. 

Eventually the full transient capabilities of the MWA 
will: 1) cover 0.25 steradian instantaneously, — approximately 
~3% of the visible sky — with full polarization cover- 
age; 2) search blindly over that area orders of mag- 
nitude deeper than previous work; 3) simultaneously 
record light curves and spectra for hundreds of sources 
and conduct routine observation of tens of thousands 
of sources per week; 4) obtain typical exposures of 
days/source, with hundreds-to-thousands of hours for 
sources in selected areas (Galactic plane, EOR fields) 
over one year; and 5) record full electric field time- 
series for a subset of observations. 

New sources detected by the MWA will have rel- 
atively poor (approximately arcminute) localization. 
This will hinder interpretation as many possible coun- 
terparts will be present at other wavelengths (e.g., [Ka^ 
plan et al. 2008). Prompt multi- wavelength follow-up 



one of which is an X-ray binary, but the other remains 
unclassified, although it varies on long timescales (~months) 
Later sear ches (Hyman et al. 2005| also see Hyman 



et al.|2009 1 found a much more rapidly variable source. 



with fiares on timescales of minutes that repeated (some- 
times) every 77min (GCRT J1745 in Fig.[8|. The posi- 
tional coincidence with the GC, along with the greatly 
increasing source density near the Centre, suggests it 
is at a distance of about 8 kpc, but this could be a 
selection effect because all the fields searched for tran- 
sients in the program were in that direction. [Bower[ 
et al. (|M05) found another transient near the GC, but 
they deduced that it was an X-ray binary with a jet 
outburst interacting with the dense medium in that 
region. 

For searches of high-latitude fields, [Bower et al.| 



(2007|) followed by [Fran et al.[ ( [2012| ) searched archival 
VLA data at cm wavelengths and identified at least 



with wide-field imaging instruments will be crucial to 
help mitigate this, as detection of contemporaneous 
variability at other wavelengths will significantly im- 
prove the localization of sources and will pave the way 
for more detailed follow-up observations such as spec- 
troscopy. 

Source confusion will limit the sensitivity of the 
MWA for source detection. However, time variability 
is a classic way to help beat the confusion limit (i.e., 
searches for pulsars with single-dish telescopes). Sub- 
tracting images of the sky taken over successive time 
periods will help beat down the confusion noise and 
allow transients to emerge, especially when correlated 
with observations at other wavelengths. Finally, exam- 
ination of polarised images that have much lower back- 
ground source levels will allow searches for bright and 
highl y polarized transients like GCRT J1745-3009 (jH^- 
|man et al.|2005i ). 
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4.4 Pulsars 

With approximately 2000 now known (Fig. [9]l , pulsars 
offer opportunities to study a range of astrophysical 
phenomena, from the microscopic to the Galactic. The 
MWA's wide coverage in frequency and large collect- 
ing area, makes it a unique instrument for low fre- 
quency detection and detailed studies of pulsars and 
fast transients. The areas of pulsar science that the 
MWA would probe at low-frequencies include pulsar 
phenomenology and pulsar searches, in general, and 
single-pulse and giant pulse studies, in particular. This 
is highlighted by the successful detection of giant pulses 
from the Crab Nebula pulsar at 200 MHz with only 
three prototype MWA antenna tiles ( Bhat et al.|2007 l. 

The MWA will be able to observe at rates faster 
than the 0.5-8 s correlator integrations using a volt- 
age capture system capable of recording, for several 
hours, correlator input data. This system will provide 
3072 X 10 kHz Nyquist sampled channels of undetected 
voltages from each of the 128 antennas of the the MWA 
that will enable high time-resolution studies of several 
types of objects, including pulsars and fast transients 
over the ~ 1000 deg^ tile beam. 

It is at MWA frequencies that pulsar emission peaks 
(Maron et al.||2000| iKondratiev et al.||2012|) and that 



interstellar propagation effects are easiest to study, but 
this wavelength range has been neglected for years be- 
cause of difficulties in processing the data and the 
increasing prevalence of radio frequency interference 
(RFI). This has left the low- frequency pulsar popula- 
tion largely unconstrained. Since low frequencies select 
against distant objects (especially for fast pulsars) the 
MWA will provide important new constraints for the 
local, low-frequency population. As a baseline, in an 
8 hour pointing at 200 MHz we should have a 8ct sen- 
sitivity limit of 1.8 mjy for long-period pulsars. When 
corrected to 1400 MHz using a spectral index of —1.8 
(Maron et al. 20001, this would be 0.05 mJy, compa- 



rable to or better than ongoing surveys like the High 
Time Resolution Universe survey (Kei th et al.|[20Tot 



or the Green Bank 350 MHz D riftscan Survey (Boyles 
et al.|2013| [Lynch et al.|2013[ ) 



Studies of individual pulsars help address a number 
of areas of astrophysics, providing unique insight into 
the as-yet poorly understood pul sar emission m echa- 
nism and astrophysical plasmas (Melrose 20031. As- 



sembling an average pulse profile from thousands of 
individual pulses can help investigate the emission fur- 
ther, probing the emission altitude, beam geometry, 
and inclination of the magnetic field from the spin 



axis (Jenet et al. 


2001 




Kondratiev et al. 


2012 





Gupta fc Gangadhara 



2003 



Hassall et al. 2012). Statis 



tics derived from average pulse profiles may be among 
the best probes of complicated phenomena like core- 
collapse supernova (e.g., [Johnston et al.|20 05). Longer 
timescale variations in the emission can directl y con- 
strain the density of the magnetospheric plasma ( Kram ;r 
et al.|2006[ ) 



Single-pulse studies should be possibl^for a sam- 



°The fractional magnitude of pulse-to-pulse variation in 
the conal profiles is expected to be close to unity. Tliere- 
fore, although desirable, it is not crucial to have the ability 



pie of pulsars in excess of 100, with the full MWA sys- 
tem at frequencies > 100 MHz ( Kondratiev & the LO- 



[FAR Pulsar Working Group[[2012p . In cases of bright 
pulsars, we aim to probe the polar emission regions at 
different heights in greater details to study single-pulse 
fluctuations simultaneously across a wide range of ra- 
dio frequencies, believed to come from different heights 
from the star surface. Mapping of the emission pat- 
terns corresponding to different heights simultaneously 
will enable tomography of the pulsar emission cone 
dMaan et al.[[20T3| [Xsgekar fc Deshpande[[200T| [Desh-[ 
pande fc Rankin|1999[ ). The wide bandwidth of MWA 
can probe a large fraction of the emission cone, and 
at larger distances from the star than higher-frequency 
telescopes, revealing much needed clues about the mech- 
anism of pulsar radio emission, which is still poorly 
understood even after 40 years of pulsar studies. 



4.4.1 Rotating Radio Transients 

While observations of pulsars with regular steady pulses 
reveal much about those sources, it is in sources where 
the pulsar emission becomes irregular that we can hope 
to gain disproportionate insight about the pulsar emis- 
sion mechanism and the relation of pulsars to other 
more exotic classes of neutron stars. In particular, the 
rotating radio transients (RRATs; [McLaughlin et ah] 
2006 and Figure |9| likely serve as links to other more 



extreme populations. The RRATS emit only sporadic 
dispersed radio bursts of between 2 and 30 ms in dura- 
tion that are separated by 4min to 3hr, but nonethe- 
less are identifiable as rotating neutron stars by their 
long-term spin-down. Initial estimates suggest that 
there could be as many RRATs as traditional pulsars 
( [McLaughlin et al.|2006l [Keane et al.|2010| ), but their 
properties do not appear the same, as they occupy dif- 
ferent parts of the P — P diagram (see Figure l9| . 

The study of RRATs and sources with and without 
radio emission may help link the canonical rotation- 
powered pulsars — neutron stars that emit primarily in 
the radio and that are powered by rotational energy — 
to other classes of neutron stars that are powered via 
other energy sources. In particular, the X-ray and 7- 
ray sources known as magnetars ( [Woods fc Thomp-[ 
json 2006), which manifest observationally as the cur- 
rentljj^ll confirmed anomalous X-ray pulsars (AXPs) 
and 9 soft gamma-ray repeaters (SGRs), are powered 
by decay of superstrong (> 10^* G) magnetic fields. In 
contrast, the soft X-ray sources known as isolated neu- 
tron stars (INS; Kaplan|2 008l are powered by residual 
heat, but may be the decayed remnants of old mag- 
netars ( [Kaplan fc van Kerkw ijk 2009). There is cur- 
rently a great deal of interest in linking radio pulsars to 
these populations. At least two magnetars have been 



seen as transient radio pulsars (Camilo et al. 2006 



2007 1, while radio searches for the INS have so far been 



to detect individual pulses from the sequences for assessing 
primary parameters of these fiuctuations, since the detec- 
tion of relatively stable fluctuation features/properties ben- 
efits from the total integration time in quite the same way 
as the average profile estimation does. 

'^http:/ /www. physics. mcgill.ca/ pul- 
sar /magnetar / main.html 
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negative ( Johnston 2003 Kaplan et al. 2003 Kon- 



dratiev et al. 2009^ . The populations may be linked 



via RRATs: the X-ray and timing characteristics of 
the RRATs offer tantalizing clues to their relation to 



both magnetars ( Lyne et al. 


2009|) and the INS ( Popov 


et al.|2006 McLaughlin et a 


.|2007 1, but as yet nothing 



Identification of RRATs relies on detecting single, 
dispersed pulses rather than steady repeating pulses 
(as in |Keane et al.|2011|[Burke-Spolaor fc Bailes|2010| ). 
They do appear as repeating transients and so many 
of the techniques required will be similar to those for 
other transients (Lorimer et al. 2007 Keane et al. 2012 1. 
We will a l so exp lore the bispectrum techniques of |Lawl 
& Bower ( 2012[ ), which enable computationally simple 
searches for dispersed pulses at the expense of signal- 
to-noise, although the small area of each antenna leads 
to a low signal-to-noise ratio on each baseline, making 
such searches ( Law et al.|2012 1 significantly less sensi- 
tivity than traditional beamforming. 



4.4.2 Interstellar Medium 

Beyond the study of pulsars themselves, variations in 
individual pulses with time and frequency can be used 
to probe the properties of the IS M ^Hassall et al.|2012 \ 
and the Galactic magnetic field ( Noutsos et al.||2008p , 
as well as the solar corona ( Ord et al.|2007 \ . Signals of 
short intrinsic duration do not survive passage through 
interplanetary, interstellar, or intergalactic plasma un- 
scathedjRTckctt 1990} [Cordes et al.| |199H [Taylor fc 



Cordcs 1993 ) and are distorted via dispersion ( [Cordes 
^Lazio„2002|), sc intillation ( [Rickett||1969J , scattering 
( Bhat et al.|2004 i, and Faraday rotation. 



Dispersion makes lower frequencies arrive later. Com- 
pared to a signal at infinite frequency, a signal of fre- 
quency u arrives with a time delay (x u^^. Dispersive 
delay can be both an advantage and a disadvantage. 
A radio counterpart to a transient event, such as a 
gamma-ray burst, observed with a high-energy satel- 
lite, would be delayed by several minutes to hours be- 
fore it arrives in the MWA band. This delay is ben- 
eficial since it allows time to prepare for the event at 



low frequen cies (e.g., Cordes fc McLaughlin|2003 Mac 
quart|2007 l. But the converse is also true. If the MWA 



provides the first detection of a transient signal, it may 
be too late for other facilities at higher frequencies to 
observe it. 

Scattering is the result of multi-path propagation 
in ionized media. While it has other efi'ects, such as 
spectral and angular broadening, we are primarily con- 
cerned with temporal broadening here. This leads to 
information on time scales shorter than the scatter- 
ing time being lost, and this time scales as cx i/"^^''^. 
Because of this steep scaling this tends to set a min- 
imum useful timescale for low-frequency observations. 
The scattering timescale becomes comparable to the 
standard MWA correlator integration time of 0.5 sec- 
onds in the Galactic plane for objects beyond ~ 5 kpc 
( Cordes & Lazio 2002). Hence, the MWA will be most 



Solar, Heliospheric, and Iono- 
spheric Science 



The MWA will make significant contributions to the 
fields of solar, heliospheric, and ionospheric science 
through a diverse set of novel radio measurements, 
as described here and in earlier papers ( Bastian 2004 



Cairns|2004|[Oberoi fc Kasper|2004| [Salah et al.|2005 
Oberoi fc Benkevitch[[20To K These measurements di- 
rectly address space weather, defined as the study of 
the impact of the Sun on the Earth and its neighbor- 
hood. Space weather effects are mostly driven by solar 
fiares and coronal mass ejections (CMEs) that pro- 
duce electromagnetic radiation (from X-rays to radio 
wavelengths), energetic particles, shock waves, changes 
in magnetic fields, and fast-moving plasma outflows 
that couple to the magnetosphere, ionosphere, neu- 
tral atmosphere, and surface of the Earth. These re- 
sponses can directly affect humanity's environment, 
technology, and society ( Scherer et al. [2(305 Warmuth 
'fc Mann 2004). The MWA will be able to observe and 



track CMEs, shocks, and other plasma structures from 
points close to their inception through their journey to 
the Earth and beyond, and it will be able to charac- 
terise their impact on the terrestrial ionosphere, mak- 
ing it a powerful instrument for space weather studies. 



5.1 Solar Science 

Even the quiet radio Sun is a time-variable, spectrally 
and morphologically complex, extended source. Emis- 
sion features are especially remarkable during periods 
of high solar activity. The active Sun produces both 
broadband (A/// ~ 1) and narrowband (A/// <C 
10~^) emissions at frequencies from below 1 MHz to 
approximately 30 GHz ( e.g., [Wild et al.|1963| 
et aL|[l998| [Cane et al ' 2002 Karlicky 2003 
2011 1. Variations in the emission can occur over time 



Bastian 
[Cairns 



scales ranging from less than a millisecond to the eleven- 
year cycle of solar activity. The complexity of so- 
lar radio emission presents a challenging radio imag- 
ing problem, but the unprecedented instantaneous dy- 
namic range, wide-bandwidth, high spectral resolu- 
tion, and full-Stokes capability of the MWA will enable 
spectroscopic imaging of the dynamic Sun. 

The physical sources of solar radio emission in the 
MWA band fall into three distinct catagories: 1) ther- 
mal emission produced by hot electrons in all layers of 
the solar atmosphere, but especially the solar corona 
(McLean fc Labrum 19851; 2) gyro-synchrotron and 



synchrotron emission with broad power-law spectra in 
frequency, created by electrons accelerated in magnetic 
loops and by shocks driven by coronal mass ejections 
fBa stian et al.|[2001 l; and 3) intense bursts of nearly- 
coherent narrow-band emissions that are produced as 
the so-called type I, II, HI, IV, and V bursts (e.g.. 



Wild et al.|1963| [McLean fc Labrum|T985l [Cane eTal] 



sensitive to scattering efi^ects for sources within this 
distance. 



2002 Cairns 20111. The presence of mulitple differ- 
ent emission mechanisms, along with propagation ef- 
fects like refraction and scattering, endow radio ob- 
servations with a rich and unique coronal diagnostics 
capability. 
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5.1.1 Electron Density and Temprature 

Large gradients in electron density and temperature 
from the solar chromosphere to the coronal regions 
have the consequence that the optical depth to radio 
waves tends to be a sensitive function of the ray path, 
which itself is determined by the refraction due to 
the density gradients. The optical depth of ray paths 
at different frequencies is dominated by contributions 
from different coronal heights. Spectral imaging by the 
MWA, therefore, should make it possible to construct a 
three-dimensional distribution of electron density and 
temperature in the lower corona for the quiescent Sun. 

5.1.2 Electron Acceleration and Radio Bursts 

Solar radio bursts are one of the primary remote signa- 
tures of electron acceleration in the corona and inner 
heliosphere. A wealth of information about bursts has 
been amassed, primarily from studies of their emis- 
sion in frequency-time dynamic spectra. Their imag- 
ing studies remain rare. The standard models for these 
bursts involve generation of radiation near the electron 
plasma frequency, Vpe, as well as near 2vpe, where: 

Upe ~ %nj'^ MHz , (4) 

for electron number density. Tie, in m~^. This mech- 
anism involves generation of electrostatic Langmuir 
waves and their conversion into ra dio emission via vari- 
ous linear and nonlinear processes (Ginzburg fc Zhe lez-| 
makov 1958, Melrose 1980; Robinson & Cairns 2000f 
At MWA frequencies, beams of en- 



Cairns 



20111 



ergetic electrons typically move away from the Sun, 
hence emission frequency is expected to decrease with 
time since the coronal electron density decreases with 
radius. 

We will focus on type II and III bursts with the 
MWA. Type II bursts are associated with suprather- 
mal electrons accelerated at magnetohydrodynamic shocks 
to speeds of 400-2000 km s~^, whereas type Ills are 
due to bursts of electrons moving at speeds of order 
one third of the speed of light, prod uced in magnetic 
reconncction sites during solar flares ([Wild et al.|1963 
■Melrose 1980; Nelson & Melrose'T985 Basti an et"aL] 
1998, , Robinson fc Cairns,2000, Cairns 2011) . 

Several major unresolved questions exist. The first 
is whether type II bursts are associated with a blast 
wave (i.e., a freely propagating shock) or the bow wave 
found upstream of a CME that is expanding into the 



heliosphere ( Cane fc Erickson 2005 Vrsnak & Cliver 



2008 [Cairns 20111. The second problem is that es 



timates of density profiles that are derived by apply- 
ing the plasma frequency relation to time-dependent 



spectra (Cairns et al. 2009 Lobzin et al. 20101, have 



tended to be much shallower than expected, as shown 
in Figure [To] The density profile is coupled to the for- 
mation mechanism of the corona, hence discrepancy 
with models casts doubt on the source of energetic 
electrons. Lastly, unresolved theoretical issues exist 
due to the relatively small levels of circular polariza- 
tion observed in type II and III bursts, despite existing 
theories predicting essentially 100% circular polariza- 
tion in the sense of the o mode (Icfthand) (Nelso n fc| 
Melrose|1985| [Suzuki fc Dulk|1985t|Cairns,2011) . 



The rapid spectroscopic imaging capability of MWA 
will address these issues by localizing burst emission. 
An initial example of this capability is shown in Fig- 
ure |11[ in which the MWA development system im- 
aged a solar burst at 15 arcminute angular resolution 
The 2 arcminute angular resolu- 



(Oberoi et al. 



20111 



tion of the full MWA will be nearly an order of mag- 
nitude better. This will help to determine the shock 
mechanism by providing radio positions to compare 
with CME observations (e.g. from STEREO), and it 
will improve density profile estimates by eliminating 
the need to infer burst location based on plasma fre- 
quency. In combination with the recent detailed mod- 
els for the source regions and dynamics spectra of type 
II burst emission, as reviewed by Cairns [ ( [2011 1, the 
MWA images and dynamic spectra should help to de- 
termine definitively the relationship between type II 
bursts, CMEs, and shocks and to provide strong tests 
of our understanding of the propagation of events be- 
tween the source and Earth. 



5.1.3 Direct Imaging of Emission from CME 
Plasma 

Though CMEs are routinely imaged in visible light 
using coronagraphs, and CME-driven shocks produce 
broadband radio emissions over several hundred MHz 
to as narrow as 20 kHz, the CME plasma itself has been 
difficult to observe directly. Meter and decimeter fre- 
quencies are regarded as the most suitable in which to 
try to image the synchrotron, gyrosynchrotron, and/or 
thermal bremsstrahlun g emission from the CME plasma 
(Ba stian fc Gary|1997 1 . Such measurements offer very 
useful and unique diagnostic capabilities, including the 
ability to constrain the CME magnetic field remotely. 
To date, however, only a few events have been confi- 



dently interpreted this way (e.g. Bastian et al. 2001 



Maia et al. [2000 1. Observing these events is partic 



ularly challenging for existing instruments due to the 
large contrasts between the inherent brightness of the 
non-thermal emission from the CME plasma (weak) 
and solar disc (strong), and other phenomenon which 
are often simultaneously present (e.g. noise storms), 
and the large angular scales of this relatively weak 
emission. With its significantly higher dynamic range 
imaging ability and higher sensitivity to extended emis- 
sion scales, owing to its numerous short baselines, the 
MWA should be the most suitable instrument yet for 
direct imaging of emissions from the CME plasma. 

5.2 Heliospheric Science 

With an increasing reliance on the robust functioning 
of satellites, humankind is now more affected by space 
weather than ever before. Reliable predictions of space 
weather are needed to protect our technological assets 
in orbit as well as humans in space. The largest of the 
space weather events are caused by CMEs, and the 
geo-effectiveness of CMEs is closely tied to the pres- 
ence of a southward directed magnetic field in the CME 
or behind the CME-driven shock. Currently the only 
operational means of obtaining this information is from 
the in-sttu observations by the near-Earth satellites. 
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Figure 10: (a) Dynamic spectrum of the type III burst of 09:26 UTC on 20 October 1999, observed 
by the Potsdam- Tremsdorf Radiospcctrograph, with normaUzed intensity (color coded) as a function of 
frequency and time ( Cairns et al. 2009, ). Horizontal bands are due to RFI. (b) Red diamonds show the 
{f,t) locations of flux maxima for each time. The green line is the best fit to ^{t) — a{t — b)^^^^ and 
implies that ne(r) ~ (r — 1.0i?s)~'^ with /3 — 2.0 ± 0.3. This unexpected result is a solar wind-like profile 
for ne(r). Fits for the leading power-law term /? = 6 of the standard density model based on coronagraph 
data (blue line) and an exponential gravitational-settling profile (black line) agree poorly with the type 
III data. 



This approach suffers from two significant disadvan- 
tages: 1) These measurements provide a lead time of 
only a few hours, consequently the advance warning 
comes so close to the event itself that it leaves lit- 
tle opportunity to take any remedial measures. 2) At 
these crucial times, interactions between the CME's 
magnetic field and plasma significantly alter the cou- 
pling between Earth's ionosphere and magnetosphere 
and the solar wind, leading to major changes in en- 
ergetic particle fiux, plasma density, plasma motion, 
magnetic and electric fields, and ionospheric activity. 
Degradation in satellite hardware and communication 
performance then occur, as well as many other space 
weather effects. 

5.2.1 Faraday Rotation 

Observation of Faraday rotation of background polar- 
ized radio sources is the only known remote sensing 
technique sensitive to the orientation and magnitude 
of the interplanetary magnetic field. Changes in FR as- 
sociated with the passage of the solar corona or CMEs 
through the line of sight to spacecraft were convinc- 
ingly detected using the telemetry signal from the Pi- 
oneer 6 and Helios satellites fStelzri ed et al.|1970 Bird 
[1982: Bird ct al._1985,) . More recently, coronal FR has 
been detected using VLA observations of radio galax- 
ies as the Sun passes in front of them ( Mancuso & 
Spangler|2000{|Spangler fc Whiting|2009[ ). The ability 



lar corona and the solar wind, including to determine 



to routinely monitor and track CME magnetic fields 
would revolutionize the field of space weather. 

Faraday rotation measurements have been applied 
to many interesting studies of the nature of the so- 



the mean magnetic field in the corona (Patzold et al 



1987 1, to measure the power spectrum of coronal waves 
and turbulence as a function of distance ( [Efimov et al.| 



T996| , to sugges t the existence of inte rmittent waves at 
single locations ([Chashei et al.|1999 l, to support theo- 



retical models of wave dissipation as a coronal heating 
mechanism (HoUweg et al. 19821, and to study the 



structure of the heliospheric current sheet (Mancuso 
fc Spangler|2000| ). 



Studies so far have been severely limited by the 
paucity of suitable linearly polarized sources and the 



narrow fields of view of existing instrumentation ( Span- 
igler fc Whiting 2009) . Observations using the WSRT 
at 340-370 MHz have lead to detection of polarized 
extra-galactic sources exhibiting typical polarized in- 
tensities of 20 mjy and readily measura ble RMs ( Wierii iga 
et al.|1993 Haverkorn et al.|2003a|c|d I . In the few ar- 



eas of the sky investigated in detail, these efforts reveal 
about one suitable polarized source every ~4 square 
degrees, as well as patches of diffuse polarized emis- 
sion ( Bernardi et al.|[2013 ), suggesting that there will 
be sufficiently numerous background sources to chart 
CMEs. 

The combination of the sensitivity of the MWA, 
its low frequency of operation, and its wide field of 
view, along with the availability of a large number of 
polarized sources with sufficient flux density, present 
the possibility of being able to measure the coronal 
and heliospheric FR simultaneously to several different 
lines of sight and out to much larger distances from the 
Sun (~ 50 solar radii) than previously accessible. No 
measurements of FR due to CME plasma are available 
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Figure 11: The top panels s how a set of radio i mages at 193.3 MHz from the MWA development system 

From left to right, the first one comes from a time close 
'10 s and occupying the entire ^30 MHz of 

-5 s and --5 MHz 



on the same intensity scale (Oberoi et al. 2011) 
to the peak of a weak non-t 



thermal emission feature lasting 
observing bandwidth, the second is during one of numerous weaker events that spanned - 
and last one shows an image of steady solar emission for the nearly quiescent Sun, although with an active 
region still visible in the northern hemisphere. The dynamic range of these images is ^2500 and exceeds 
that of earlier images by about an order of magnitude. The bottom panels show spatially localized spectra 
of emission across the solar disc at the same times as the top panels. The spectra are shown for every 
third pixel, the pixel size is 100" x 100". The spectra span the entire observing bandwidth, binned into 24 
frequency bands spaced 1.28 MHz apart and averaged over 0.8 MHz. The y-axis ranges in arbitrary units 
are (left to right) 50-4620, 100-6000 and 100-1800. Celestial north is on top, and the red circles mark 
the size of the optical solar disc. 
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yet at large solar offsets, hence we rely on simu l ations 
to build quantiative expectations. Jensen et al. ( |2010 1 
estimate the RM from CMEs to be ^0.05 rad m~^ 
(12° at 150 MHz) at about 0.4 AU. |Liu etal] ( |2007 l 
also arrive at a similar estimate and also show that the 
FR time series from instruments like the MWA during 
the CME passage may be used to infer the helicity and 
size of the CME. 

For multiple reasons, these are rather challeng- 
ing measurements to make ( Oberoi fc Lonsdale|2013 1. 
Ionospheric RM, is 1-10 rad m~ , more than an order 
of magnitude larger than the expected heliospheric or 
CME RM signature. This implies that we will need to 
achieve ionospheric RM calibration accuracies of a few 
percent over large fields of view in order to study helio- 
spheric RM. The quiet solar flux density is 10'*-10^ Jy, 
hence detecting 20 mjy sources requires polarimetric 
imaging dynamic range of about 10^ , in the vicinity of 
a strong and variable Sun. 

5.2.2 Interplanetary Scintillation 

As radio waves from a distant compact radio source 
propagate through the turbulent and inhomogeneous 
solar wind plasma, the incident plane wavefronts pick 
up phase corrugations corresponding to density inho- 
mogeneities. These phase corrugations develop into 
intensity fluctuations by the time the wavefronts ar- 
rive at Earth-based telescope. Due to the motion of 
the solar wind, the interference pattern produced by 
these intensity fluctuations is swept past the telescope 
giving rise to observed scintillations. 

The evolution and dynamics of solar wind in the in- 
ner heliosphere dominates space weather effects. How- 
ever, there are few techniques available to study the 
solar wind in this vast region. Historically, interplan- 
etary scintillation (IPS) was the first remote sensing 
technique used to investigate the solar wind in this 
region ( [Hewish et al. 19641. More recently, routine 
availability of observations from heliospheric imagers 
(e.g., STEREO) that image very large parts of the 
heliosphere using white light scattered from from the 
solar wind electrons, have emerged as very useful tools 
for investigations in the inner heliosphere. These tech- 
niques, along with local in-situ sampling provided by 
most satellite based instrumentation, complement each 
other well since they cover different spatial scales and 
have lead to a more complete characterisation of the 
solar wind than wo uld have been possible with either 
of the datasets (e.g. Jackson et al.|20"07 Dorrian et al. 
[2008l l. 

Several new instruments are coming online to study 
IPS. Apart f rom the MWA, the new facilities inc lude 
the LOFAR ( |Bisi et al.|2011| [Fallows et al.||2012[) and 
Mexico Array Telescope ( [Mejia-Ambriz et al.|2010l ) . In 
addition, the Ooty Radio Telescope in India (Manoha- 
ran, P. K., private communication) and Solar Wind 



Imaging Facility (SWIFT) in Japan ( Tokumaru et al. 
2011 1 are both undegoing considerable upgrades which 



observable, the normalized co-spectrum obtained from 
simultaneous observations at frequencies differing by a 
factor of ~2 (Scott ct al. 1983 , Zhang_ 2007 ). Under 
the assumption of weak scattering, the interplanetary 
medium can be regarded as a collection of thin screens 
placed along the line-of-sight, and the IPS observable 
can be modeled as a weighted sum of the contributions 
from each of the thin screens. The physical parmeters 
to which IPS is sensitive are: 1) the velocity of the solar 
wind; 2) the level of electron density fluctuations, Sue, 
3) the spectral index of the Sue power spectrum, 4) 
the inner scale for the power law governing 5ne; and 
5) the axial ratio of Sue- 

Using an inversion process, it is possible to extract 



the global 3D distribution of IPM parameters ( 


Jack- 


son et al.| 1998 |Kojima et al.| 1998| |Asai et al.| 


1998 


Oberoi|2000| 


B si et al.|2008||2010||Dorrian et al. 


2010 


Xiong et al. 


20111. The most important bottleneck 



currently constraining the fldelity and resolution of in 
version of IPS data, and hence its science output, are 
the sparseness (both in space and time) and the non- 
uniformity of the sampling of the heliosphere afforded 
by the present generation instrumentation (Oberoi & 
Benkevitch||2010| ). The MWA will represent a signifi- 
cant advance in addressing the bottleneck of limited 
heliospheric sampling when its voltage capture sys- 
tem is completed. The MWA will have a sensitiv- 
ity of ~70 mJy (1 second, 30.72 MHz), comparable 
to that of the SWIFT transit facility currently under 
commissioning in Japan. The MWA will boost the 
number of independent IPS measurements by up to a 
factor of ~30 over what is possible with most present 
day systems. The southern hemisphere location of the 
MWA also complements the heliospheric coverage pro- 
vided by the existing facilities, all of which happen 
to be in the northern hemisphere. Combining these 
data with other MWA observations like solar imaging 
and FR observations, and the observations from other 
sources like STEREO, SDO and HINODE, wiU pro- 
vide a very comprehensive dataset for charaterizing 
the inner heliosphere, monitoring and tracking CME 
evolution from close to the Sun to our terrestrial neigh- 
borhood. Further, the software framework being de- 
veloped for handling the MWA IPS data could pro- 
vide a common data repository, interface, and analysis 
tools for IPS data from observatories across the world 
( Oberoi fc Benkevitch|2010 1. 



will substantially improve their scientific capabilities. 

For the MWA, the primary IPS observable will be 
the power spectrum of observed intensity fiuctuations, 
although we a will also employ a less commonly used 



5.3 Ionospheric Science 

Propagation through the dynamic and structured mag- 
netoionic ionospheric plasma leaves an imprint on the 
incident low radio frequency radiation from astronom- 
ical sources. In order to arrive at the true proper- 
ties of the incident radiation, the ionospheric distor- 
tion must be identified and removed. The MWA can 
tackle this long standing challenge by measuring the 
refractive shift in the apparent position of dozens of 
bright and relatively compact sources at a cadence of 8 
seconds. The ionospheric calibration information thus 
comprises of a time series of vectors showing the re- 
fractive shift of calibrator sources from their fiducial 
positions. To first order, the magnitude of the shift 
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is directly proportional to the ionospheric gradient in 
the total electron content (TEC) towards the fiducial 
sky position of the calibrator source. This information 
for all the observed calibrator sources leads to a series 
of differential TEC maps and an accuracy of ~0.01 
TECU (1 TECU= lO^*^ electrons m'^). These maps, 
produced routinely as a part of the MWA calibration 
will form the primary data products for ionospheric 
science. The MWA beam will sample a region about 
175 km X 175 km at a height of ~ 300 km and the ex- 
pected calibrator source density corresponds to spatial 
scales of ~ 10 km at this ionospheric height. 

The high-precision and dense sampling of the iono- 
sphere provided by the MWA offers new opportunities 
for ionospheric science. The three areas of ionospheric 
research that the MWA will make significant contribu- 
tions are in: 1) the measuring of waves in TEC associ- 
ated with travelling ionospheric disturbances (TIDs); 
2) the development of storm enhanced density (SED) 
and sharp TEC gradients associated with space weather 
events; and 3) the tracking of the onset and evolu- 
tion of scintillation associated with the development 
of ionospheric irregularities. 



to influence the development of SED near the MWA 



f Co ster et al.|2007| [Foster fc C oster 2007| [Foster fc Er- 
|ickson||2013p . SEP associated gradients have a prac- 
tical consequence, as they cause large errors in GPS 
range measurements, and significantly degrade the high 
precision positioning information provided by GPS sys- 
tems. The MWA's longitude (~ 116°E) is very close 
to the longitude being extensively instrumented in the 
Asian sector as part of the Chinese Meridian project 
(~ 113°), and the conjugate longitude of ~ 64°W is 
within the field of view of the US sponsored chain of 
incoherent scatter radars, in particular Millstone Hill. 
The combination of simultaneous measurements from 
these distributed sensors thus provides a unique way 
to determine mesoscale response to geomagnetic dis- 
turbances in the coupled geospace system. In partic- 
ular, it will be possible to investigate whether there 
is a seasonal and universal time dependence of geo- 
magnetic storm response by comparing American and 
Asian/ Australian sector response to solar wind and 
CME forcing. Such effects have been predicted the- 



oretically ( Foster fc Erickson 2013) but have not yet 
been confirmed by observations. 



5.3.1 Traveling Ionospheric Disturbances 



TIDs are manifestations of middle-scale ionospheric 
irregularities arising as response to acoustic gravity 
waves. TIDs are frequently observed at high and mid- 



dle latitudes (Brist owet al.,1 996 NicoUs & Heinselman 
2007| [Oliver et al.||1994| |1997| |pjuth ct al. 2004|, and 



their activity and amplitudes vary depending on lati 
tude, longitude, loca l time, season, and solar cycle (e.g . 

Hernandez-Pajares et ah] 20061. 
They can be generated by auroral sources at the high 
latitude thermosphere, the passage of the solar termi- 
nator (Galushko et al. 20031, and storms, hurricanes 



and tornados in the troposphere. Detailed information 
about the sources, energetics, and scale sizes of energy 
coupling from the lower atmosphere to the upper at- 
mosphere, along with resulting ionospheric effects, is 
critical to developing a deep understanding of overall 
upper atmospheric energy balance. Measurements of 
TIDs at the VLA have been made since the early 1990s 



([Jacobson fc Erickson|1993 Jacobson et al.|1995[) and 
recent observations have been reported by |Helmboldt| 



et al 



( 2012|?[ ). The fine grid over which the MWA cal- 
ibration solutions will provide differential TEC mea- 
surements will provide detailed measurements of the 
propagation of TIDs over the array and TEC gradi- 
ents or irregularities. 

5.3.2 Storm Enhanced Density 

The MWA is well-sited to contribute to coordinate 
international investigations of the SED phenomenon 
(Foster 19931. This process is responsible for some 



of the largest space weather effects at mid-latitudes 
( Coster & Foster 2007 1 since it produces severe, tempo- 



rally dependent electron density gradients. These gra- 
dients will be obvious in the differential TEC measure- 
ments made by the MWA. The geographic/geomagnetic 
pole offset in the southern hemisphere is anticipated 



5.3.3 Ionospheric Scintillation 

Ionospheric scintillation can be defined as the rapid 
fiuctuation of radio-frequency signal phase and/or am- 
plitude, generated as a signal passes through the iono- 
sphere. MWA observations using voltage capture ca- 
pabilities will be uniquely suited for the exploration of 
ionospheric irregularities and scintillation. A voltage 
capture system that can record raw antenna voltages 
from the full array for 24 hours will be added to the 
MWA by the end of 2013. This system should provide 
the ability to create array-beams in post-processing 
that use the entire bandwidth and collecting area of 
the MWA, and that can be pointed independently any- 
where within the MWA field-of-view. Such a system 
will provide scintillation measurements at multiple po- 
sitions across the visible ionospheric region. 

MWA array-beam observations will be augmented 
by data from specially designed GPS receivers at and 
near the MWA site that can measure L-band scintil- 
lation statistics. GPS TEC data from these receivers, 
when combined with data from a large number of newly 
deployed GPS receivers on the Australian continent, 
will provide important contextual information about 
the regional distribution of TEC around the MWA 
site. Other satellite based remote sensing platforms, 
including COSMIC, DMSP, and JASON, will also be 
available during overfiights to provide further informa- 
tion about background ionospheric plasma and electro- 
dynamic conditions in the ionosphere. A similar ap- 
proach has been followed previously in the Combined 
Radio Interferometry and COSMIC Experiment in To- 
mography (CRICKE T) campaign utiliz ing the VLA 
(Dyniond ct al. 200"8| [Coker et al.[[2009| ). This mix of 
data types allows for ionospheric specification from the 
global scale down to the regional scale. 
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6 Conclusion 

The MWA is a unique low-frequency instrument with 
high survey efficiency and arcminutc angular resolu- 
tion. It will provide deep integrations on EoR tar- 
get regions that could enable the first detection of 
the redshifted 21 cm signal and detailed characterisa- 
tion of diffuse Galactic emisison and faint extragalactic 
source populations. The planned full-Stokes southern 
hemisphere sky survey will include the Galactic Cen- 
tre, Galactic plane, and the Large and Small Magel- 
lanic Clouds, making the MWA particurly well-suited 
to identifying the missing population of SNRs, and 
thereby constraining the total energy budget of the 
Galaxy. 

Long integrations in the EoR fields and regular ob- 
serations along the Galactic plane will enable sensitive 
searches for pulsars and transient events on timescales 
from seconds to days. Blind surveys for both rare and 
faint events are planned, as well as targeted obsreva- 
tions of known transient sources, including X-ray bi- 
naries, neutron stars, pulsars, and low-mass stellar ob- 
jects. 

The MWA will open new windows into space weather 
and solar physics and provide new measurements of so- 
lar bursts throughout their journey from the surface of 
the Sun to the Earth and on to the outer Solar System. 
High-dynamic range spectroscopic imaging will show 
the origin of bursts in the solar corona, while monitor- 
ing of interplanetary scintillation will probe the inter- 
stellar wind. As events reach the Earth, the MWA's 
calibration system will yield detailed measurements of 
ionospheric disturbances traveling over the telescope. 

The telescope is complementary to other new low- 
frequency radio telescopes, falling between LOFAR and 
PAPER in both sensitivity and angular resolution. It 
is the only low-frequency SKA precursor and is lo- 
cated on the very radio-quiet site planned for the SKA 
low-band telescope. The combination of ASKAP and 
MWA at the MRO will yield a comprehensive program 
of precursor science for the SKA over the comining 
years. The MWA posseses a versatile design with ex- 
cess infrastructure capability that will allow the tele- 
scope and its users to embrace new scientific and tech- 
nological opportunities over the coming decade, hope- 
fully providing an enduring and fruitful future for the 
instrument. 
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